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INTRODUCTION 


Several  programs  have  been  developed  for  plotting  contours 
on  different  types  of  plotting  devices.  Most  of  the  programs  require 
ordered  arrays  or  other  restrictions  on  the  data  being  processed.  The 
data  for  this  program  have  the  restriction  that  the  elevation  above 
the  x-y  plane  not  be  the  same  as  surrounding  points;  i.e.,  several 
points  in  an  area  should  not  be  the  same  elevation. 

The  program  either  develops  or  reads,  as  input  data,  triangular 
patterns  from  nodal  values  input.  From  these  triangular  elements, 
each  side  is  investigated  to  determine  if  one  of  the  desired  contours 
fall  in  this  element.  When  all  elements  have  been  investigated, 
the  contours  are  ordered  and  spline- filled  for  a smooth  curve  output. 

Note  that  the  density  or  sparsity  of  nodal  points  determines  the 
accuracy  of  the  contour  plots. 


II.  PROGRAM  DESCRIPTION 

The  program  is  developed  in  three  overlays.  Each  overlay 
uses  a dynamic  core  allocation  scheme. 

The  (0,0)  overlay  determines  user  preferences  as  to  plotter  type, 
number  of  points  used  in  spline- fitting  spacing  of  contour  labeling, 
number  of  contours  desired,  and  their  valves.  This  overlay  remains  in 
core  during  the  entire  process. 

The  (1,0)  overlay  reads  the  input  nodal  values  and  their  corre- 
sponding elevation  and,  if  not  specified  by  the  user,  generates  the 
triangular  elements  from  the  nodal  data. 

Overlay  (2,0)  locates  the  desired  contours  as  specified  in  the 
main  overlay  and  orders  these  contours  for  plotting.  This  routine 
produces  an  output  of  nodal  points,  elements,  and  contours  in  addition 
to  tape  7 which  is  used  in  the  next  overlay  for  plotting. 

Overlay  (3,0)  reads  the  data  from  tape  7,  spline  fits  the  data  to 
four  times  as  many  points  as  input,  then  plots  the  data  on  the  desired 
device. 

When  the  triangular  elements  have  been  developed,  it  is  only 
necessary  to  implement  the  proper  overlay  for  additional  or  different 
contour  output.  The  output  on  the  Houston  Instruments  is  plotted  on  a 
20-  X 20-inch  plot. 
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PROGRAM  USAGE 


Input  data  may  be  of  one  or  two  forms.  If  the  user  does 
not  supply  triangular  elements  as  input,  the  input  format  is  in  free 
formal  as  follows: 

CARD  1 JMAX,  (Number  of  nodes  in  plot) 

CARD  2 thru  JMAX 

X coordinate,  Y coordinate,  Z elevation 
If  a triangular  pattern  is  included,  the  format  is: 

CARD  1 JMAX,  NELM,  (Number  of  triangular  elements) 

CARD  2 thru  JMAX 

X coordinate,  Y coordinate,  Z elevation 
CARD  3 JMAX  +2  thru  NELM 

Node  1 Node  2 Node  3,  of  each  element 

The  input  is  on  Tape  1. 

When  the  program  begins  execution,  the  following  information  will 
appear  on  the  screen: 

FORHAT 

ONLV  X,Y  AND  Z COORDINATES  ARE  INPUT 
COORDINATES  AND  ELENENTS  ARE  INPUT 
REPLOTTING  OF  GENERATED  CONTOURS 
SPACING  OF  CONTOUR  NUNBERING 
SMOOTHING  FUNCTION  FOR  SPLINE  FIT 
SET  TO  I FOR  HOUSTON  INSTRUMENT  PLOTS 
SET  TO  I FOR  TEKTRONIX  OUTPUT 
NUMBER  OF  CONTOURS  TO  BE  PLOTTED  20MAX 
TO  CONTINUE 
TO  END 


NHflE 

IFLO 


INUM 

IPTS 

IHIP 

ITEK 

NCON 

CONT 

END 


UALUE 

1 

1 

2 

3 

20 

0 

0 

0 

0 

0 


U 


Regardless  of  the  input  used,  IFLO  must  be  set  to  1,  default  until  • 

a successful  run  has  been  completed.  After  a successful  run  (i.e.,  | 

tape  7 has  been  generated) , IFLO  set  to  2 may  be  used  to  generate  addi-  | 

tional  contours.  If  IFLO  is  minus,  no  contour  nodes  will  be  plotted 
(Appendix  B,  Figures  B-2  and  B-3'). 

The  following  parameters  may  be  changed  with  each  execution: 

INUM  Factor  which  controls  the  spacing  of  the  numbering  of  the 
contours  - default  20. 

IPTS  Smoothing  parameter  for  the  spline- fit  routine  - default  4. 

IHIP  Set  to  1 for  Houston  Instruments'  plots. 

ITEK  Set  to  1 for  Tektronix  plots. 

NOON  Set  to  the  number  of  contours  desired. 

CONT  To  continue  after  all  other  parameters  have  been  set. 

END  To  end  program  execution,  the  input  format  is  A4,  12. 

When  the  desired  parameters  have  been  entered  and  CONT  is  entered, 
the  program  will  ask  for  NCON  contours  to  be  entered.  These  are 
floating  point  numbers  in  free  format.  Tape  6 contains  all  of  the 
generated  output. 

IV.  PROGRAM  CHANGES 

One  subroutine  is  written  in  compass.  This  routine  returns 
the  amount  of  core  available  to  the  program.  This  may  be  replaced  by 
the  following  subroutine: 

SUBROUTINE  KOREFL  (A,  KORE , KFL) 

DIMENSION  A(l) 

C 

C SET  KOPE  equal  to  or  greater  than  35  times  the  number 

of  nodes 

C 

KORE  = XXX 

RETl’RN 

END 
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V.  HARDWARE  REQUIREMENTS 

The  following  hardware  requirements  are  necessary: 

a)  Houston  Instruments'  plotter  or  equivalent. 

b)  Tektronix  401X  terminal  or  equivalent. 

VI.  SOFTWARE  REQUIREMENTS 

The  following  software  requirements  are  necessary: 

a)  Houston  Instruments'  plotting  package  or  equivalent. 

b)  Tektronix  Terminal  Control  System  (TCS). 

Two  sample  problems  are  presented.  Appendix  A presents  a list  of 
the  input  data  for  each  sample  problem.  Appendix  B shows  the  output 
obtained  on  the  Houston  Instruments'  plotter  for  each  example.  Appen- 
dix C presents  a list  of  the  output  obtained  for  the  first  example. 
Appendix  D is  the  FORTRAN  listing  of  the  program. 
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Appendix  A. 

SAMPLE  PROBLEMS  INPUT  LISTINGS 
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INPUT  DATA  SET  I 

105  0 

fl.5  .5  3.4 

17,4  .4  4.7 
il.«  1.?  4.1 

9.6  ?.6  3.0 

16,3  3.1  4.6 
17,  3.7  4,7 
7.9  4.0  2.7 

12.2  4.9  3.9 

12.6  7.4  4.2 

7.3  8.2  3,5 

2.5  8.5  3.4 

5.1  11.4  4.3 

1.8  11.8  4.8 

16.7  12.  5.7 
18.  13.6  5,8 

15.5  13.4  5.0 

6.2  16.8  5.2 

3.8  17.9  5.5 

13.5  18.  5.9 

2.2  19.4  4.9 

o.O  20.  5.8 

17.5  19.5  5.9 

6.7  21.1  6.1 

2.0  23.2  4,6 

3.2  22.9  4.9 

4.5  23.4  5.0 

9.7  24,4  4,4 

24.4  1.4  5. 

29.?  1.7  4.1 

32.4  0.1  3.3 

29.0  3.3  4.0 

38.3  1.8  2.2 

41.8  ?.9  ?.B 

21.4  6.5  4.9 

31.3  5.0  3.3 

36.1  4.7  2.7 

41.1  5.3  3.3 

43.2  5,5  3.8 

49.4  3.8  2.8 

33.4  7,6  3.3 

22.2  10. 4.9 

27.1  10.3  4,2 

31.8  9.8  3.5 

43.3  8,5  4.1 

31.1  12.2  4.3 

30.8  13,8  4.5 

31.8  14, C 4,7 

29.8  15.5  4.8 

23.1  16. C 5.1 

94.9  12.7  4.5 

41.4  12,8  5.0 

43,3  15.4  5.8 

34.0  17.5  5.8 
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37. q IB. 3 6.4 

39.1  18,0  6.3 

33.7  19.8  6.7 

39.  c 6.7 

30.?  ?1.4  6.8 
37, c 22.0  6.8 
42. c 22.3  6.8 
?3.D  ?2.3  6.0 
??,-»  26.1  6.8 

24.0  26.3  7.0 
42. r 25.4  6.7 

48.5  25.1  5.4 

10.3  25.6  4.8 

8.7  27,0  4.5 

4.3  27.2  4.0 

2.1  29,7  3.3 
14.?  29.0  5.1 

19.7  28.4  6.5 

33.0  ’8.5  7,9 

29.6  30.5  7.7 
25.^  30.7  8.1 
48, p 30.8  4,9 

44.1  31.3  5.3 

46.0  32.9  4.7 

45.3  33.1  4.7 

39.8  32.5  5.8 
13. P 34.5  5.0 
10. P 35.8  5.1 
19,?  36,0  6.5 

20.3  36.9  6.4 
27. P 36.3  6.5 

40.  Q 36.6  4.4 

47.0  36.4  4.0 

1.1  39.4  3.9 
13.^  39.3  5.6 

33.0  38.5  5,5 

18.8  42.1  5.4 

35.3  4^.6  4.9 

38.3  41.7  4.5 

43.7  41.1  3,5 

28.8  42,2  5.2 

32.7  42.8  4.8 

15.1  44.6  4.9 

14.0  45.9  4.6 

18.7  47.4  4.2 

15.7  48.3  4.1 
18. (C  46.0  4,6 

24.0  48.5  3.9 

77.4  48.5  3.6 

38.4  45.6  3.6 

41.0  47.0  2.9 

47.0  47.5  2.7 
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INPUT  DATA  SET  2 

157  0 

0,  0.  4. 


0. 

10. 

5. 

0. 

50. 

6. 

0. 

80. 

0. 

40. 

8. 

0. 

>=:0. 

9. 

0. 

60. 

10. 

0. 

70. 

11. 

0. 

°0. 

12. 

0. 

OO. 

13. 

0. 

100. 

14. 

0. 

120. 

15. 

0. 

130. 

16. 

0. 

150. 

17. 

0. 

180. 

18. 

0. 

200. 

19. 

-3. 

ii . 

5. 

-14 

. 37 

. 6. 

-17 

. 41 

, 7. 

-21 

. 65 

. 8. 

-27 

. 84 

. 9. 

-30 

. 93 

. 10 

-33.  102,  11. 
-36.  112.  13, 
-3«5.  121,  15. 
-42.  131.  17. 
-46.  140.  18. 
-53.  158.  !<?. 
-6.  9.  5. 

-24,  33.  6. 
-30.  41.  7. 
-36,  48.  8, 
■42.  57.  10, 
-48.  65.  12. 
■54.  73.  14. 
■60.  «1.  16. 
•66,  89,  18, 
■78.  105.  19. 
■8.  7.  5. 

■16.  13.  6. 
■32.  26.  7. 
•40.  32.  8. 
•48,  T7.  10. 
•56.  44.  11. 
•64.  51.  13. 
•72,  57.  14. 
•80.  64.  16, 
■88.  69.  17, 
■96.  76.  18. 
■10^..  81.  19. 


-65. 

24. 

12. 

-74. 

27. 

14. 

-83. 

30. 

16. 

-‘?3. 

34. 

17. 

-102. 

, 38, 

. 18 

-IIP. 

. 41, 

. 19 

-10. 

-1. 

5. 

-?o. 

-3. 

6 

-39. 

-4. 

7. 

-49. 

-5. 

8. 

-6.  <3. 

-7,  !0. 
-7fl.  11. 

-«7.  -9.  12, 
-9ft.  -1,1,  13. 
-10ft.  -11.  IS. 
-lift.  -12.  IS. 
-12S.  -13.  19. 
-<9.  -4.  S. 

-?7.  -12.  6. 
-45.  -20.  7. 
-54.  -24.  8, 
-ft3.  -24.  9. 
-72.  -32.  10. 
-81.  -3ft.  11. 
-OO.  -38.  14. 
-100.  -42.  16. 
-109.  -44.  19. 
-6.  -«.  5, 


-3ft. 

-48. 

6. 

-4?. 

-56. 

7. 

-4P. 

-64. 

8. 

-84. 

-72. 

10. 

-60. 

-80. 

11. 

-66. 

-88. 

12. 

-72. 

-9ft. 

14. 

-78. 

-104. 

15. 

-84. 

-112. 

16. 

-90. 

-120. 

17. 

-9ft. 

-128. 

18. 

-10?. 

-13ft 

. 19 

0.  -10.  5. 

0.  -30.  ft. 

0.  -60.  7. 

0.  -90.  4. 

0.  -no.  9. 

0.  -no.  10. 

0.  -170.  11. 


r ^ 


0.  -IflO.  12. 
0.  -190.  13. 
<*.  -9.  5. 
n.  -18.  7. 

12.  -27.  8. 

21  -46.  9. 

S8,  -126.  10. 
f>7,  -144.  11. 

79.  -171.  12. 
10.  -3.  S. 

19.  -6.  7. 

28.  -9.  8. 

37.  -12.  9. 
67.  -21.  10. 
117.  -36.  11. 
136.  -42.  12. 
165.  -51.  13. 
183.  -57.  14. 
10.  .7  5. 

20.  1.5  6. 

30.  2.2  8. 

50.  3.7  9. 

70.  5.2  10. 
100.  7.5  11. 
110.  8.4  12. 
120.  9.1  13. 
130.  10.  17. 
160.  12.1  18. 
180.  13.5  19. 
9.  5.  5. 

26.  18.  7. 

35.  20.  8. 

44.  25.  9. 

53.  30.  10. 
62.  38.  11. 
70.  40.  12. 

80.  45.  13. 
97.  55.  15. 
105.  60.  17. 
139.  80.  18. 
173.  100.  19. 
5.  P.  5. 


13. 

27. 

6. 

18. 

36. 

7. 

23. 

45. 

8. 

27. 

54. 

10. 

32. 

63. 

11. 

40. 

81. 

12. 

44. 

90. 

13. 

49. 

99. 

15. 

83. 

108, 

, 16 

58. 

117, 

. 17 

66  . 

135. 

, 18 

80. 

162. 

, 19 

Appendix  B. 
CONTOURS 


NUMPFO  OF  CONTOI'OS  TO  PF  PLOTTFO 
MUMOFO  OF  NOOFS 
MljMnFD  OF  FIFMFMTS 


F 

’ OF 
! P? 


^OOF 

x-'“O0Po 

v-COOPO 

i 

.PSOOPOOEtOl 

.5PP0PP0r*00 

? 

.174COOOE*OP 

,4eopooo‘"*oo 

7 
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.1200PP0b*O1 

<4 
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.17POnoOE*PP 
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7 
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a 
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n 
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IP 
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• 1140000‘^*02 

n 

.1POOOOOE*01 
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,1F70000E»0? 

.1200000F*02 

.lPO0O00E*0? 

,13F0000F*02 

IF 

.1SFOOOOE*0? 

.1340000F*02 

17 

.F2POOOOE*01 

.16BOPOOF*0? 

IP 

.OPOOOOOE^OI 
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lo 
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.1POOPOOP*0? 

?n 
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PI 

.PPOOOOOE*01 

.2000000F»0? 

?p 

.175PPP0E*0P 

.1950000'^*02 

?3 
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.211OPOO‘^*02 

2u 

.2000000F*01 

.2320000'^  *02 

?= 
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?o 

.2P20P0OE*P2 

.1700000f*01 

30 
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31 

.2900000E*02 

.33POOOOF*01 

3? 

.3P30P00E«02 

.1P00000‘^*01 

30 

,41POoooE»0? 

.250COOOB  *01 

3^ 

.P140O00E+0P 

.F500000P*01 

3F 

.3130n00E*P2 

.5000000F*01 

3F 

.3610000E*02 

.47000000*01 

37 

.4HOOOOE«02 

.5300000'-*01 

3P 
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.55P0P00F*01 

3° 
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.3POOPOOF*P1 
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.7FO00C0F*01 
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.1220000'-*02 

PF 
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.13R0O00‘^*02 

47 
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4P 

.2PP0P00F*O2 

.15500000*02 

4P 

.231"O00F»i»2 

.1FP0000B*02 

SO 

,3400pnof *P2 

.12700000*02 

51 

.414000PE*P2 

.12800000*02 

S? 

.4330P00E»02 

.15400000*02 

F3 

,3400O00E*02 

.1 7500000*02 

54 

.37POooOF‘02 

. lP30»>0Po  *02 

V&l  LF 


.taOOPOPF*'' 
,470C000F*P 
.il00n00F*n 
.OOOCOOnE*0 
,4^00P00E*p 
.■iTOf  cooE*n 
.?7PonppE*o 
.7qOOCOflE*fl 
,4?PCPOOE*o 
.OFP0000E*P 
•0400POOE»n 
.^3000nOE*0 
.4«nopooE*« 
,p70ono9E+n 
,';soonooF»o 
,sqooPOOE*p 
.'^?P000PE*P 
.‘:EP0P00E*n 
.SPOOPOOE^P 
.^.PCPPOOF*0 
.FBOonooE^p 
.sqocnooE*p 

.‘>10CP00E*P 
.iFOOPPOE»P 
,t,Q0000OF*0 
.cOO0n00E*0 
.4AOOOOOE  »0 
.FOOOOOOE*0 
.4iocpooE»n 

.3300P00E*fl 
.POPOOOOE»n 
.??OOPOOE*P 
.?POOOOOE»0 
.4P00P00E*0 
.33O0POOE*P 
.?700000E*0 
. 33PeP00E*P 
.3POOPOOE*0 
.?POOPOOE*0 
.33OOP0OE*0 
.aPPP000E*0 
•4?COPOOE*0 
•OcoOPOOE^O 
.6100P00E»0 
.6300P00E*P 
.4F00000E*0 
,^.700P00E*P 
.4p00P00E»n 
.CIOOOOOE^P 
,450PPP0E*P 
.<=0P000PE»P 
,PBPPPP0F*P 

,sPorpooE»n 

,<1*P0P0PF  *0 


^0 


I 


5*" 

.3910O00E *0? 

. IPOOOOO'^02 

.6300000E+0 1 

. 3370''00E*0? 

.19RPP006^0? 

.6700P00E+P1 

^7 

,39cpooOE>0? 

,?0S00006»02 

.6700000E+01 

5P 

.30P0O0OE40? 

.2l40000r*02 

.6800n00E+01 

SO 

,07'^OOOOE*OP 

,??PO000F+0? 

.6POOOOOE+01 

ftr. 

.4250000E*0? 

.??30000'"^0? 

.6800POPE+01 

61 

.?3P''OOOE^O? 

.?230P00r-4.0? 

.6000000E+01 

h? 

.2330O00Ef0? 

.2610000c-*0? 

.68POOOOE+P1 

60 

,?4OOp00E»0? 

,2630000<-402 

.7000000E+01 

6^^ 

.4?O000OE*0? 

.?540000'**02 

.6700POOE+01 

6S 

.4flS0O0OE+O? 

. 25100000+0? 

.S400P00F+01 

6S 

,in30O00E*0? 

.256000CF+02 

.4800P00E+01 

^7 

.8700000E^01 

.2700000r-+02 

.4500000E+01 

f-  P 

.430a000E^Ol 

.27?COOO'^*02 

.4POOPOOE+01 

60 

,?1POOOOE*01 

.297PO00F+0? 

.3300POOE+P1 

70 

.14?OOOOE^O? 

.2900000'^  + 02 

,S100000E»01 

71 

.1970000F>0? 

.2840000F+0? 

.6500P0OE+01 

7? 

,3300000E»C‘2 

.2850000P+02 

.70POOOOE+01 

70 

.?O6O000E^0? 

,3050n00t-*0? 

.77O0POOE+P1 

7u 

.?S60000E*08 

,3070000F^O? 

.PIOCOOOE+Pl 

7<^ 

.4P?0O00E*O? 

.3080000F+02 

.4900000E+P1 

76 

,441 OOOOE^C? 

.3l30000r»0? 

.S300P0PE+01 

70 

.46O0''00E*0? 

.3290P00r+02 

.4700000E+01 

7P 

.4SO0op0E^0? 

.3310P00>^^02 

.47P0P00E+01 

70 

.39SOPO0E»O9 

.32sppoO'-+0? 

.SPOOPOOE+01 

«n 

.13POp00E^0? 

,3450000F^02 

.SOPCPOOE+Pl 

PI 

,IOPO«OOE*0? 

.358PP00F+0? 

.FlOOPOOE+01 

p? 

.19?0O00E*O? 

.3600000- ♦O? 

.6'^OOPOOE  + Pl 

P o 

.2030000E^O? 

,36OP000r+n? 

.6400P00E+01 

84 

.?7ROO00E*09 

.363OP00F+02 

.6'=^POPPOF+01 

PS 

,4090000F40? 

•3660000O.0? 

.44POPOPF+01 

8^ 

,47O0P00E»0P 

,3640000!'  +02 

.4C00P00E+P1 

87 

. 1 lOOOOOE^Ol 

.3940000' ♦0? 

.OO0C900E+P1 

pp 

.13?0000E^0? 

.3930000'^*02 

.S600000F+P1 

oo 

.3300000F^0? 

,385OO0O(-  + 0? 

.5SOOPOOE+01 

00 

.IRROOPOE^O? 

.S?1PP00!-*02 

.S4P0P0OE+01 

91 

.OS3O000E*02 

.4060000'+02 

.4Q0PP0PF+ni 

07 

.36OOP00E*07 

.417O00C--»02 

.4«!:POPOOE+C1 

93 

,437r)op0E*0? 

.4110OOOF+0? 

.ofooPOOE+01 

0/* 

,7PPPO00E+0? 

,42?PP00p  +02 

.s?opnpOF+oi 

QC 

.3?3PO00F*o? 

.4280000'^^02 

.4P0CP00E+01 

0<^ 

.l'=^ino00E*0? 

,4460000''  + nP 

.49000POE+01 

97 

. 140PPOOE*0? 

,4590OP0''*0? 

.46PP000E+ni 

OP 

,1P3o«(,0E*0? 

,474PPPO!-+02 

.42POPOOE+81 

OO 

. ISOOPPOE^O? 

.4830PP0O+02 

.41O0P00F+P1 

inr 

, lP60000E^n? 

.460POPOF»O2 

.460PPP0E+01 

101 

.?400nOOE^O? 

.4850P00F+P2 

. 390»'00PE+Ol 

loo 

.3?40«OOE^O? 

,485PP0Pr  *02 

,3600000F+ni 

100 

,3P4PopOP*OP 

.4S60P0Pt'»0? 

. 3600PPOE+P1 

104 

.41S0000E»0? 

.470000PF+02 

,?900POPE+P1 

ino 

.4700000E^O? 

.475POOOr  +02 

.2700P00E+01 

21 


9R 

68 

?6 

25 

99 

26 

27 

23 

100 

<S9 

68 

24 

101 

17 

12 

18 

102 

12 

16 

10 

103 

21 

17 

23 

104 

9<=; 

102 

92 

105 

PI 

95 

92 

106 

«5 

91 

92 

107 

76 

85 

78 

105 

96 

97 

99 

109 

100 

96 

98 

no 

101 

100 

98 

ill 

34 

6 

2 

112 

42 

31 

35 

113 

31 

34 

28 

114 

9 

10 

8 

115 

6 

9 

8 

115 

45 

43 

50 

117 

43 

42 

40 

U5 

76 

75 

77 

119 

75 

64 

65 

120 

64 

60 

65 

121 

60 

59 

57 

122 

57 

54 

55 

123 

55 

50 

51 

124 

26 

25 

23 

125 

25 

68 

24 

126 

12 

1« 

13 

127 

18 

17 

23 

125 

10 

12 

11 

129 

102 

Q7 

103 

130 

92 

R9 

93 

131 

55 

78 

86 

132 

78 

76 

77 

133 

99 

96 

98 

134 

6 

? 

5 

1 3=^ 

2 

34 

28 

136 

31 

35 

29 

137 

35 

42 

40 

13« 

26 

31 

29 

139 

10 

A 

7 

140 

8 

6 

5 

141 

43 

50 

40 

142 

75 

77 

86 

143 

60 

65 

5? 

144 

57 

60 

55 

145 

50 

44 

1 46 

51 

55 

52 

147 

25 

23 

20 

24 


148 

24 

25 

20 

149 

18 

13 

20 

150 

13 

1? 

11 

1=^1 

23 

1« 

20 

152 

11 

10 

7 

153 

92 

103 

93 

154 

103 

10? 

104 

155 

85 

93 

86 

156 

78 

86 

77 

157 

2 

5 

3 

158 

35 

29 

30 

159 

40 

35 

36 

1^0 

8 

7 

4 

161 

5 

8 

3 

16? 

50 

40 

36 

163 

65 

52 

39 

164 

52 

60 

55 

165 

51 

44 

52 

166 

44 

50 

37 

167 

103 

93 

104 

168 

93 

36 

105 

16Q 

30 

35 

36 

170 

7 

4 

1 

171 

4 

8 

3 

172 

36 

c;o 

37 

173 

52 

38 

44 

174 

37 

48 

38 

17>^ 

93 

104 

105 

176 

36 

30 

32 

177 

4 

1 

3 

178 

37 

36 

32 

179 

39 

44 

38 

180 

38 

37 

33 

181 

32 

37 

33 

1®2 

38 

3Q 

33 

25 


rONinUR  NUMREP 
ELEMENT 


1 CONTOUR  VALUR  = .3000000E+01 

LINE 


139  .89750000F+ 

15?  .76750000E* 

15A  .40?OOOOOE* 
150  .33700000F* 

160  .R9750000R* 

16?  ,34750000R^ 

163  .AR993333F* 

167  .A1R66667P* 

16P  .47P07692F* 

169  .33700000F* 

170  .95983017F* 

17?  .35900000F* 

173  .4fl993333F» 

175  ,41R66667F* 

176  ,34?50000R* 

178  .38600000F* 

170  .48461538F* 

180  .41520000F* 

181  .40336364F* 

18?  .48160000R* 


01  .4??50000F*01 

01  .55750000E*01 

02  .47214?86E*02 

02  .48500000E+01 

01  .4??50000E*01 

02  .61500000E*01 

02  ,45733833E*01 

02  .46016667E+0? 

02  .44938462E*02 

02  .48800000E*01 

01  ,26013O86E*01 

02  ,60333333E*01 

02  ,4S733833E*01 

02  ,46016667E*0? 

02  .24000000E*01 

02  .50000000E+01 

02  .45230769E401 

02  ,36200000E*01 

02  ,48454545E*01 

02  .41400000E+01 


,76750000E*01 

,55857143E*01 

.41057143E+02 

.34750000E+0? 

,95983017E*01 

.35O00000E+02 

,49330769E*02 

.41057143E+02 

.45762500E+0? 

,34?50000E*0? 

.8157l4?9E*01 

,38600000E*0? 

.48461538E*02 

.45762500E+0? 

.34O09091E*0? 

,40836364E*0? 

.48160000E+0? 

.42080000E*0? 

.41520O00E*0? 

,42080000E«0? 


.557O0000F*01 

.59285714E+01 

.46800000E*02 

.61=00000F*01 

,26013986F*01 

,60333333E*01 

.543R461SE*01 

,46800000E»02 

,45100000E*02 

,24000000E*01 

.25000000E*01 

,50000000F»01 

.45230769E*01 

.45100000E*02 

,56363636E*00 

.43454S4SF+01 

.41400000E*01 

,31000000F*01 

.36?00000E*01 

.31000000E*01 


26 


CONTOUR  NUMREP 
ELEMENT 

X 


RB  .19O03333E* 
6?  .23653333B* 
74  .54R33333r* 
7P  .19083333E* 
79  .18117ft47E* 
BO  .323666f>7E4 

BP  .23653333F4 
87  ,?4369?31B4 
90  .110B5714E4 
9?  .31362'500F4 
100  .42907436F* 
10?  .90083333F* 
104  .3P366667F4 
110  .23228571F4 
11?  .29001314F4 
114  .110B5714E* 
IK-  .12333333F* 
Ilf  .31362S00F4 
117  .28442857F* 
12P  .59^50000F♦ 
120  .34133333F4 
130  .42144444E* 
13f  .?9001314f* 
137  .20O33333F* 

140  .12B00000F4 

141  .33350000F* 
1*50  .42333333F* 
1S3  .37466fif7F4 
155  .42144444F* 
15P  .29462500B4 
161  ,12785714F* 
16?  .34275000F* 
163  ,46960000F* 
166  .37483333b* 
16P  .47P96643F* 
171  .12000000F* 
17?  .35233333F* 

173  .46960000B* 

174  .43025000F* 
177  ,11328571F* 
179  .43769231b* 


? CONTOUR  VALUF 


Y 


01  .3Q100000E*0? 

02  .4B073333E*0? 
01  .3?072??2E*0? 
01  .39100000E*02 

01  .39394118E*0? 

02  .46600000E*0? 
02  .4P073333E*0? 
02  .48015385E*02 
02  .76?85714E*01 
02  .11300000E*02 
01  .27201428E*0? 

01  .9?P33333E*01 

02  .46600000E*0? 
02  .4P142857E*0? 
02  .33009709E*01 
02  .76285714E*01 
02  .57333333E*'U 
02  .11300000E*0? 
02  .10157143E*02 

01  .10?OOOOOE*0? 

02  .45477778E*02 
02  .3B600000E*0? 
O'.  .33009709F*01 
02  .91222?22E*01 
02  .47E00000E*01 
02  .11?50000E*0? 

01  .10433333E*0? 

02  .43P66667E*0? 
02  .3B600000E*02 
02  .21125000E*01 
02  .4f428S71E*01 
02  .10575000E*0? 
02  .84400000E*01 
02  .96166667E*01 
02  .36403757E*0? 
02  .30500000E*01 
02  .1047777BE*02 
02  .84400000E*01 
02  .P1000000E*01 
02  .11000000E*01 
02  .81384fl5E*01 


= .4000000E*01 

LINE 

X 


.18il7B47E 
.?4369?31E 
.42987436E 
.54P33333E 
.29428571E 
.31FOOOOOE 
.23?28S71E 
.31B00000E 
.90083333E 
.28442857E 
.20461R38E 
.59P50000E 
.34133333E 
.22100000E 
.2P033333E 
.12333333F 
.12POOOOOE 
.333500C0E 
.28500000E 
.42333T33E 
.37466667E 
.40000000E 
.29462500E 
.28E00O00E 
.12785714E 
.34?75000E 
.22000000E 
.40000000E 
.47B96643E 
.29600000E 
.12000000E 
.35233733E 
,489P4f 15E 
.43025000E 
.47P99723E 
.1 IfOOOOOE 
.374P3T33F 
.43769331E 
.43?66667F 
.1 IBOOOAOE 
.43266667E 


Y 


♦01  .3939411BF*02 
*0?  .4P0153B5E*0? 
*01  .?7?014?fiE*02 
*01  .32072???F*0? 
♦01  .403?R571E^02 
♦0?  .469?5000E^O? 
♦02  .48142857E^02 
♦02  .469?5000E^02 
♦01  .92P33333E^01 
♦02  .10157143F^0? 
♦01  .26200000E^0? 
♦01  .10?OOOOOF^02 
♦02  .4547777rE^02 
♦0?  .48133333E^0? 
♦02  .9122????F^01 
♦02  .57333333E^01 
♦02  .47500000E^01 
♦02  .112‘=000OE^02 
♦0?  ,97000000E^01 

♦ 01  .10433'»T3F^0? 
♦0?  .43P66667E^02 
♦02  .41400000E^02 
♦02  .21 125000E^01 
♦O?  .97000000E^01 
♦0?  .464?B571E^01 
♦0?  .10‘=75000E^02 
♦01  .9914?PF7E^01 
♦0?  .41400000E^02 
♦0?  .364037R7E^0? 
♦0?  .lEOOOOOOE^Ol 
♦0?  .30500000E^01 
♦0?  .1047777flE^02 
♦0?  .1363076QF^02 
♦0?  .PlOOOOOOF^Ol 

♦ 0?  .364034UF^02 
♦02  .132727?7F^01 
♦0?  .96166667E^01 
♦0?  .813P461BE^01 
♦0?  .75000000E^01 
♦0?  .13?727?7E^01 
♦0?  .75000000F^01 


CONTOUR  NUMRER 


T CONTOUR  VALUE 


EooonooE*oi 


element 


17 

X 

.766000007*0? 

13 

.1?9000C07*0? 

19 

.766000007*0? 

71 

.796400007*07 

7? 

.179000007*07 

73 

. 1 14090917*0? 

7° 

.?11333337*0? 

31 

.793333337*0? 

34 

.30?109?6E*0? 

34 

. 13?633337*0? 

36 

.106816187*0? 

4= 

.717333337*07 

44 

.774500007*07 

46 

.743833337*07 

5? 

.304400007*0? 

59 

.!0450000=*0? 

97 

.117700007*07 

96 

.999144477*01 

99 

.309900007*07 

40 

. 349909097*0? 

41 

.148?89717*0? 

4? 

.701844677*0? 

43 

.119I?9007*0? 

4 = 

. 119647067*0? 

71 

.374000007*0? 

7? 

.349538467*0? 

73 

,440990006*01 

74 

.107090917*0? 

79 

.99959996=*01 

77 

.940000007*01 

76 

.000166677.01 

79 

.891941 187*01 

60 

. 109900007*0? 

61 

.379000009*0? 

6? 

.369909097.0? 

»3 

.403000009.0? 

64 

.140769717.0? 

69 

.196400007.0? 

04 

.701666677.0? 

67 

.180619389.07 

6“ 

.706I79009.0? 

90 

. 1 19647069.07 

91 

.147866679.07 

94 

.471750009.0? 

91 

. 3=6894749.0? 

90 

,449990009.01 

99 

.450437977.01 

lOl 

.999555567.01 

lO? 

.96=000007.01 

I01 

.410333337.0? 

113 

.741850797.07 

116 

.4117=0007.07 

119 

.470959967.07 

173 

. 360666677.07 

174 

.449353737.01 

114 

.434166677.01 

13? 

.447000007.0? 

139 

.743850757.07 

I 30 

.7440??997.07 

14  = 

.41  3935067.07 

141 

. 149166677 .01 

149 

.137147869.01 

151 

.746666679.01 

143 

.449766679.0? 

145 

.414010=59.0? 

111 

,449766677.0? 

LINT 

< 


.?9o<,rtoooE«0? 
.Ilf-ORORIE  *0? 
.?ST33T3TF*0? 
.30?10R?»'E«0? 
.13’«33'»3E*0;? 
, JOf>Hl«lflE»0? 
,??fS00#0E*0? 
.?A3«3333E*0? 
.30E400COE»0? 
.1045P000E«0? 

.n??0000E«02 

.?1‘^i?';ooe*o? 

.?3S444AAE*P? 
,?3';44444E  ^0? 
.3?400e00F*0? 
.10?090<51E»0? 
.R4000000E  »01 
,«9?94nfiE«01 
.3?>:OOrOOE*0? 
.40300000E»0? 
.lSP4ftn00E*0? 
,?R061=;3RE*0? 
.?OP1?>=99E»0? 
. 147RE4fi7E  *0? 
.34RR3o44E*0? 
.15f,R9i.74E»0? 
,45043?97E*OI 
.103164E7E«n? 
.96COPOOOE  *01 
.RE4117EEE*01 
.IP3I6467E*0? 
.134flCOO0E*0? 
.?9?50CC0E*0? 
.34o16447E*0? 
.34olf,Nf,7E*0? 
,43o33733E«0? 
. I 34R0000E  *0? 
.lft700OO0E*0? 
. Ifl-’OOPOOE  *07 
.79?<;O0OOE*O? 
.lE9inO00F*0? 
.174?5n00E*0? 
•IRolOIOOE^O? 
.47oSSRR6E*07 
.36nAfi4ft7E  *0? 
.aaoTRtpTE  *01 
.RN4ll74EE*«l 
.43<,164f)7E  *01 
.l?4?5Pfl0F*0? 
.44700P00E  *0? 
.?<,/.0?7R9F*0? 
.4S=OOflO0E*O? 
.4R7EOOOflE*0? 
.4n93904E*0? 
.TAolft4(,7F*01 
,?3714?RftF*01 
.4S=00P00E  *0? 
.74  7«R7R3E*07 
.?440?NESE*0? 
.41401ORSF  *0? 
.PRTROOPOF  *01 
,?44444ft7E*rtl 
.7R7SOOOOF  *01 
.4P43P44?F  *07 
.431000(107  *0? 
,431000007*07 


.140333337*0? 
.77(1464477*07 
.14433333E*07 
.14090000E*0? 
.?7846467E*0? 
.?4490909E*0? 
.110777777*07 
.1R"33333E*07 
.IS957637E*07 
.7R444647E*0? 
.?43lfilB7E*0? 
.110777777*0? 
.13350000E*0? 
.1S433333E*0? 
. 199000007*0? 
.34333333E*07 
.716400007*0? 
.361000007*07 
.47<500000F*0? 
.374343647*07 
.41Ra?«57E*07 
.436064677*0? 
.106679007*07 
.10??39?9E*0? 
.149949457*07 
.149461947*0? 
.?3401999E*0? 
.390181677*0? 
.196000007*0? 
.7?914?64E*07 
.341000007*0? 
.193357047*0? 
.419000007*0? 
.419714797*07 
.376361647*0? 
.346470977*0? 
.4364?«97F*0? 
.441000007*0? 
.436064677*0? 
.431697317*0? 
.716790007*01 
,lfl??39?9F*0? 
,90C33-,5if  .01 
.309790007*07 
.141734647*0? 
.73401699E.O? 
.734006337*0? 
.196000007*0? 
.177790007*0? 
.330644477*0? 
.14?537317*01 
.309750007*0? 
. 309000007*0? 
.14177???7*0? 
.731979I?E*0? 
.191914677*0? 
.3?>000007*0? 
.147531317*01 
.140094997*01 
.177999007*0? 
.??7900O07*0? 
.139476977*0? 
.191900007*0? 
.173046477*0? 
.1P79741?F*0? 
. 1 7304447F  *0? 


Y 


.160900007*0? 
.744909097*0? 
.196333337*07 
.I999?43?7*0? 
.?86666477*07 
.?4?16167E*0? 
.133900007*07 
.156333337*0? 
.199000007*0? 
.343333337*0? 
.716400007*07 
.106479007*0? 
.153664477*0? 
.153664677*0? 
,14994949E*0? 
,350161677*0? 
.??914764E*0? 
.393357947*07 
.419i14?97*0? 
.346476977*0? 
.441000007*07 
.431497317*07 
,716790007*01 
.969333337*01 
.145441947*0? 
. 141736647*0? 
.734006337*0? 
.354611117*0? 

.1??79oOOE*0? 
.?3?35?94E*0? 
.354611 117*07 
.43?60000F*0? 
.479679007*0? 
.4C?90000E*0? 
.407900007*07 
.330466677*0? 
.437600007*0? 
.440900007*0? 
.440900007*0? 

.479679007*0? 

.619000007*01 
.114=00007*0? 
.619000007*01 
. 309000007*0? 
.1417???77*0? 
.?339791?7*0? 
.?3?19?047*0? 
. 151914677*0? 
. 114900007*0? 
.3??000007*0? 
.140094997*01 
.3?100000F*n? 
.?94600007*P? 
.177999007*0? 
.??7=000n7*0? 
.1394?6=7E*0? 
.3?10000CF*0? 
,1396334)7*01 
.140014667*01 
.1?79741?7*0? 
,19941467F*0? 
. 191=00007*0? 
. 19=41 6477*07 
.?10?30777*0? 
.1?1=?Q417*0? 

.l?19?94ir*0? 


CONTOUR  NUMRER 


t*  CONTOUR  VALUE 


6000000E»01 


element 


Y 


LINE 


Y 


3 .17R66667E* 

5 ,?37624?5E* 

7 .17R6iSftft7E* 
10  .237624?5E» 

12  .23R029R9E* 

14  .17735714F* 

lA  .17400000E4 
19  .23799534E* 

21  .30040000E* 

25  .30400000E4 

27  .17400000E4 

2P  .l‘Tlfll82E* 
3?  .3?  «80952E* 

34  .32263158E* 

40  .30400000E4 

41  .32227273E* 

4?  ,2ei84615E* 

43  .15866667P* 

44  .19700000'^* 

50  .38R80952E* 

S3  .33933333E4 
6P  .40055556E* 
7?  .35300000E+ 

77  .68764706E* 

94  .43050000E* 

97  .37268421E4 

99  ,#iP764706E* 

103  .66444444F* 

119  .44411111E.4 

120  .45928571E* 

123  .3P400000E4 

124  .64083333E* 

127  .66444444E* 

143  .45R28S71E4 

14f  .39630769E4 

147  ,f.4083333E4 

151  .62166667E* 

164  .43140000E* 


02  ,?0983333E*02 

02  ,22?83300F*0? 

02  .20983333E*0? 

02  , 222833005*02 

02  .222966025*02 

02  ,286142P6F*02 

02  ,35500000E*0? 

02  .222958035*02 

02  .19040000E*02 

02  .374000005*02 

02  .3SS00000E*02 

02  .387727275*02 

02  .321190485*02 

02  .182157895*02 

02  .374000005*02 

02  .366818185*02 

02  .385692315*02 

02  .378333335*02 

02  .389800005*02 

02  .321190485*02 

02  .180111115*02 

02  .309222225*02 

02  .177666675*02 

01  .212941185*02 

02  .283500005*02 

02  .171210535*02 

01  .212941185*02 

01  .206222225*02 

02  .275000005*02 

02  .239000005*02 

02  .171168675*02 

01  .212500005*02 

01  .206222225*02 

02  .239000005*02 

02  .168000005*02 

01  .212500005*02 

01  .205666675*02 

02  .167800005*02 


.180333335*02 

.180333335*02 

.177357145*02 

.237995345*02 

.300400005*02 

.177333335*02 

.177333335*02 

.238029q95*02 

.322631585*02 

.322272735*02 

.15866667^*02 

.197000005*02 

.384578955*02 

.338333335*02 

,281846155*02 

.384578955*02 

.231333335*02 

.190181825*02 

.231333335*02 

.400555565*02 

,353000005*02 

.430500005*02 

,372684215*02 

.746666675*01 

.444111115*02 

,384000005*02 

.650000005*01 

.746666675*01 

.455000005*02 

.455000005*02 

.396307695*02 

.650000005*01 

,621666675*01 

.431400005*02 

.416200005*0? 

,632500005*01 

.632500005*01 

.416200005*02 


.202333335*02 

.202333335*02 

.286142865*02 

.222958035*02 

.190400005*02 

.304333335*02 

.304333335*02 

.222966025*02 

.182157805*02 

.366818185*02 

.378333335*02 

.389800005*02 

,322894745*02 

.180111115*02 

,385692315*02 

.322894745*02 

,386666675*02 

,387727275*02 

.386666675*02 

.309222225*02 

.177666675*02 

.283=00005*02 

.171210=35*02 

.207333335*02 

.275000005*02 

.171166675*02 

.213090915*02 

.207333335*02 

.252384625*02 

.252384625*02 

.168000005*02 

.213090915*02 

.205666675*02 

.167800005*02 

,164400005*02 

.209583335*02 

,209=83335*02 

.164400005*02 


9 


rONTOUP  NUMREP 

S CONTOUR  value 

= ,7000000E*01 

ELEMENT 

LTMF 

X 

Y 

X 

1 

.24R90932F*02 

,26299302E»0? 

,24P92730E*02 

2 

,24R90932F*02 

,26299302E*02 

,24R99?31E*02 

4 

,24R92730F*02 

,26302936E*02 

•21543750E402 

6 

.24P9923lF*02 

,26297202E»0? 

.249ia4R5E*02 

9 

,21543750F»02 

,29118750E*02 

.21?00000E*02 

ii 

.24918485F*02 

,26282910E*02 

,30066667E*02 

IF 

,28550000F*02 

,33883333E*02 

.27il25«0E*02 

17 

,21200000F*02 

.34343750E402 

.22i70588E*02 

20 

,30066667F»02 

,23422222E*02 

,30709091E*02 

25 

.28550000F*02 

,33883333E*02 

,30681R18E*02 

26 

.27112500F*02 

,34550000E*0? 

,22i705fl«F402 

32 

,35785714F*02 

.302142R6E402 

,33?47368E*02 

33 

,30709091r*02 

,22690909E*02 

.33F250C0E+02 

41 

,30681818F^02 

.33045455F*0? 

.33?47368F*02 

50 

,35785714F*02 

.30214286E*02 

,39750000E*02 

5i 

,33525000F*02 

,21975O00E»02 

,36681R1PE*02 

69 

,39750000F*02 

.26175000F*02 

.366P1P18E*02 

Y 


,?630?Q16.E*0? 
.?6?9720?E*02 
.291187S0E*02 
,26?f»2910E*02 
,343A37‘;'>F*02 
.234222?2E*02 
,34550000E*02 
.3A71176SE«02 
,22f9090QE*02 
,330454SSE»02 
,3471 176SE*02 
,31236«4?E*02 
,2197S000E*02 
,3123ft«4?E*02 
,2617SOOOE*02 
.23iaiPlPE*02 
.231fllPlPE*02 


30- 


1 


CONTOUR  NUMBER 

6 CONTOUR  VALUE  = .ROOOOOOE+01 

ELEMENT 

X 

LINE 

Y X 

Y 

4 

.?5536364r+02 

.30300000E+02 

.252312505*02 

.305562*^05*02 

n 

.?5'=:36364P  + 02 

.303000005*02 

.26600000E*02 

.306=;00005*02 

0 

.252312«=;0r*02 

.305562505*02 

.25200000E*02 

.310312=05*02 

15 

,26ft00000F>02 

.30650000E*02 

.257375005*02 

.310500005*02 

17 

.25200000F+02 

.310312505*02 

.25288235E*0? 

.310647065*02 

.?5737500E*02 

.310500005*02 

.252882355*02 

.310647065*02 



3 

COMTOUP  X 

Y 

,3000F^01 

,S5«=;714E>01 

.S9?a571F*01 

.lOOOF^Ol 

.767‘^OOOF^Ol 

.5S7«:OOOF^01 

,3000F^01 

,B97S000E^0l 

,42?'^000F*01 

.3000F^01 

,959B30?F^01 

.2601 399F^01 

,3000F^01 

,R1';7143F^01 

.250O000E»01 

-.3000F^01 

,^0?OOOOE^O? 

,4721429E^02 

-.3000F^01 

.A1057UFf0? 

.4680000E^02 

-.3000F^01 

,41«»S6F'7E^O? 

,4601667E>0? 

-.3000F»01 

,^57ft?5CE*0? 

.451f000F^02 

-,3000F*01 

,47?07f>OF^0? 

,449-»846F^02 

,3000F^01 

,^933077E^0? 

,F43o462E»01 

,3000F*01 

,489Q333F»0? 

,457'<333E401 

.3000F4-01 

,452'>077E^01 

,3000F^01 

,4818000E^O? 

,414'>000E^01 

,3000P*01 

,4a08000E*0? 

.BlOftOnOFtOl 

.30rtOF*01 

,41'^POOOE^O? 

,3623000F>01 

. 3000F+01 

,4033638F»0? 

,434q4FSE^01 

,3000F*01 

.38BOOOOF*0? 

.5000000F*ni 

.3000F*01 

,38Q0000F*0? 

.603t333F*01 

,3000F*01 

,347S000E*0? 

.615C000F^01 

.3000F*0! 

,3370000E*0? 

,485n000F*01 

. 3000F*01 

,34?=;000E^0? 

.2400000F^01 

.3000F4.01 

.1400909F*0? 

.‘^63^364E^00 

! 

32 

H 

f 

' 

■ 

CO^'TOUP 


X 


Y 


.^►OOOP^Ol 
,4000F^0l 
,4000F»0l 
,AOOO>^*Ol 
,4000F*01 
.4000F*01 
-,^OOOF*01 
-,4000F*01 
-,4000F*01 
-,4000F*01 
-,4000F*01 
-.4000F^01 
-.4000F^01 
-,4000F*01 
-.40n0F^01 
-,4000F*01 
-.4000F^01 
,4000F*01 
,4000F^01 
,4000F^01 
,4000F^01 
,4000E^01 
,4000F*01 
,4000F*01 
.4000F^01 
.4000F»01 
,4000F^01 
,4000F*01 
-,4000F*01 
-,4000F>01 
-,4000F*01 
-,40nOF^Ol 
-,4000F^01 
-,4000F>01 
-,4000F>01 
-.4000E»01 
-,4000F^01 
-.4000F*01 
-.4000F+01 
-.4000F  *01 

-.4000E*rtl 

-,4000F*01 

-.4000F»01 

-,4000F»01 


,?94?«S7F^0l 

.l8n'^f>SF^0l 

,190M31'<F^01 

,5443313F^01 

,4?9B744E*01 

,?04^1‘;4F*01 

,?2lO<iOOF*0? 

,23??8S7F^0? 

.?3^>R333F^0? 

.2436923F^0? 

,31S0000E^0? 

,32386ft7P*0? 

.3413333F*0? 

.3746667F^02 

.40nOCOOF*92 

,4?14444E^0? 

,47flQ6f>4E^02 

.2200000E*0l 

.4233333E^01 

.5925000E»0l 

.900ft333E*01 

.1108571F^9? 

.1233373F-»0? 

.12fi?009E^0? 

.)27rtS71E^0? 

.I20nOOOF*0? 

.nsooooE^o? 

,1 13?857F^02 
,2980000E^O? 
,294ft?S0E^0? 
.2900131E*0? 
,?8n3333F^0? 
,20«=nOOOF*O? 
,2044286E*O? 
,3138?F0E*0? 
,3335000E^0? 
,3427';0(1E>0? 
,3523333F^0? 
.3748333F+02 
.430?‘=:00E*0? 
.43?8687F*n2 
,4374Q^3E*0? 
,46P8000F>0? 
,489H487F^02 


,403?8'^7F»0? 
,393Q412F*02 
.39l''OPnF^O? 
,320''2?2F*02 
,?72-'143F^0? 
.?6?''000E-»02 
,481->333E*0? 
.48l62fl6F»02 
,4P07333E»02 
.4P01538E>02 
.469-5S00F*0? 
.466^000F^n? 
,4547778E*02 
.4388667E+02 
,414pOOOE»02 
,3P6r>000E^02 
.364'^376F^02 
,9qU?P(SF^01 
,in47333F*02 
.102o000E^02 
.9?8'5333F»01 
,762n=;7ir>01 
.S73-'333F*f)l 
.475''080r*01 
,464'>PS7F»01 
,305''000E^01 
,1327273F^01 
.llOoOOOF^Ol 
,150'>OOOF*01 
.21 1^500F*0l 
.330'>Q71F*01 
.912P??2F»01 
.970'’Cif>0E*01 
.101'^714F*0? 
.int'OOOF^OP 

.n?'^oooF»o? 

,1057500F^02 
,1'>477'^PF*C2 
,Q6i A887fc"*ni 
.PlOOOOOP+Ol 
.7F0'^0n0P  + Pl 
.81  3‘^4^2P*91 
.P44''000'^*0  1 
.1 38"077F*C? 


33 


COMTOUR 


X 


Y 


.RflOOP 

.ROOOF 

.=OOOF 

-.500OF 

-.SOOOF 

-,=OOOF 

-.SOOOF 

-.FOOOF 

-.SOOOF 

-.ROOOF 

-.5000F 

-.FOOOF 

-.FOOOF 

-,FOOOF 

-.FfOOF 

-.FOOOF 

-.FOOOF 

-.SOOOF 

-.SOOOF 

-.FOOOF 

-.ROOOF 

-.SOOOF 

-.FOOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.=OOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.=OOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.FOOOF 

-.=OOOF 

-.SOOOF 

-.SOOOF 

-.=OOOF 

-.SOOOF 

-.=OOOF 

-.SOOOF 

-.SOOOF 

-.■"OOOF 

-.FOOOF 

-.=OOOF 

-.FOOOF 

-.SOOOF 

-.''OOOF 

-.FOOOF 

-.FOOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.FOOOF 

-.SOOOF 

-.SOOOF 

-.SOOOF 

-.FOOOF 


01 

.?i,R0?00F‘0? 

. 14009S0F.01 

01 

.?R1HS07F‘0? 

.I42F373F.01 

01 

.?<,1RB1SF«0? 

.139a33SE*01 

01 

,i,fl?S0OOF»O7 

.39FO000F.0? 

01 

,i,7«SR5ftF»0? 

.305r000F*03 

01 

.b717S00F*07 

.3093SOOF*0? 

01 

.45S0000E»0? 

.3?1''OOOE*0? 

01 

.4O70000F*0? 

.3?2O000F*0? 

01 

.43013riF*0? 

.330fO67E*0? 

01 

.4030000F*0? 

.3<.e<.?O0F*0? 

01 

.36SQ0RIf *0? 

.376  7016E*0? 

01 

.34RIB67F*07 

.402SOOOF*02 

01 

.33SOOOOF*0? 

.41571<,3F.02 

01 

.7OSSO00F*0? 

.42S''OOOF.02 

01 

.39?SOOOF«0? 

.<,39b7SOF*02 

01 

.?fl0S154E»0? 

.<.31''Q33F*02 

01 

.?01«SS7F*0? 

.43flS607F.O? 

01 

.lS7n000F»0? 

.<.<.0'"000F*02 

01 

.15«4000F*0? 

.441f>000F*02 

01 

. 14a?PS7E*0P 

.<.30<,3BoF*O2 

01 

.UOaOOOF^O? 

.432f0O0F*02 

01 

.«97R41?F*01 

.3933S?9E*02 

01 

.R9R16f7F*01 

.361'‘OO0F.O? 

01 

. 1 031FA7F.((5 

.3S<.f11  1F*02 

01 

.10?0R09E*0? 

.3S01 «lflF*0? 

01 

.104SOOOF»0? 

.3<.37333F*0? 

01 

.1 3?S337F*0? 

.300O007F.02 

01 

.13ROOOOE»0? 

.378oo07F*03 

01 

.1160R0OF»0? 

.3440091F*02 

01 

. 10««1R?E»0? 

.2421 «1BF*02 

01 

. 1 1??OOOF*0? 

.2IP'0O0E*02 

01 

.9400000F*01 

.22S’ 439F.0? 

01 

.«6B1 1 7BF*01 

.?32’S39F*02 

01 

,4S04330E*01 

.214O0«3F*02 

01 

.<,49Qq00E*01 

.2340190F*02 

01 

.01 

.233O7S1F.0? 

01 

.3401067F.01 

.337'<000F.0? 

01 

.7S7RoOOF*CI 

. 195<.167F.03 

01 

.24FS6fr7F«01 

.191=000F*O? 

01 

.337lo?nF.01 

.n5<.2SSF*02 

01 

.43<.16ft7F*01 

.1S1S107F.03 

01 

.S9SF5FSF»01 

. 1S6O000F.02 

01 

.96S00r0F»01 

.1227SO0F.O3 

01 

.1?4?SOOF*P? 

.114a000F*02 

01 

.139FB71F.0? 

.10233S3F.02 

01 

. 147B667F*0? 

.R0S7333F.O1 

01 

.16910O0F*O? 

.oi9in00F*0l 

01 

.?0B1?S0F*0? 

.7187‘^flOF.Ol 

01 

.31S1?S0F*0? 

.I0P*'2S0F.02 

01 

.?173333F*07 

. 1 10’23?F*02 

01 

.?20=OOOF*03 

. 1 33"000F*O2 

01 

. 23S444BF  *03 

. 1S3F007F*03 

01 

.?43«333F*0? 

. ISOYl'nF.O? 

01 

.?S31333F*0? 

. lS0-n33F*O3 

01 

.?RflOOOOF*0? 

. 10S->333F.02 

01 

.?9HOOOOF*03 

.lo0''0O0F*03 

01 

.30310S1F*0? 

.159'=3f,3F.02 

01 

. 30F4000F*07 

.1S9"OOOF*02 

01 

.3?40000F*0? 

.I<,9s<,FSF.02 

01 

.34SS3«SF*0? 

.14S<,015F.O? 

01 

.3SSB9<,7F*0? 

. 141 73F8F.02 

01 

.300SF67F*03 

.1417322E*02 

01 

.ono.isiF.o? 

,127oQ90F*02 

01 

.4l4010SE*0? 

.12797MF*02 

01 

.4330000F*P3 

.121'204F*02 

01 

,440?(,S7F*03 

. I23''OF7F*02 

01 

.<,00  104, be  *03 

.31fl330flF*0P 

m 
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COK'TOUP  X 


,f.000F*01 

,6000P>01 

.FOOOF^Ol 

.^OOOP^Ol 

-,6000F*01 

-,6000F>01 

-.^OOOF^Ol 

-.FOOOF^Ol 

-,F000E*01 

-,60ft0F^01 

-,F000F>01 

-.ftOOOF^Ol 

-.ftOOOF^Ol 

-.FOOOF^Ol 

-.6000F-*-01 

-,6000F*01 

-.6000P^01 

-.IS000F4.01 

-.eoooF-^oi 

-,frOO0F*Ol 
-,F)000F>^01 
-.FOOOF^Ol 
-.6000F+01 
-,6000F*01 
-.FOOOE^Ol 
-,6000E*01 
-,6000E>01 
-,F000E*01 
-.F000F*01 
-,F000P>01 
-,FOOOE*01 
-,6000F>01 
,F000E*01 
.FOOOF^Ol 
,6000E*01 
,6000E*01 
,6000E^01 
,f OOOE*01 
,6000F*01 
,6000F*01 


Y 


. 17R(S6'S7E^O? 
.IBOmiE^O? 
,?376?41E*0? 
,?37qqS3E^O? 
,17B6667E^0? 
,177.3S71E*0? 
.1773333E*0? 
,1740000E*0? 
,15fi6667E^0? 
,19riPlflE-»0? 
.1970000E*0? 
.?313333F^0? 
.?81B46?F^02 
,304ftOOOE*O? 
,3???7?7E*9? 
,384S78qE*0? 
,38flROqFF4.0? 
,400F5‘=^E^O? 
.43O'^OOOE^0? 
.444) lllE^O? 
.45FOOOOE*0? 
.45q?P57E>n? 
,4314000E*0? 
,416?OPOE^O? 
,39P3077E*0? 
.3840000F*0? 
,37?6P4?E^0? 
,3530000E*0? 
.33q3337E*0? 
,323P316E*0? 
,3O040O0E^0? 
,23P0300E^0? 
.P87P471E^01 
,7468667E*Ol 
, 6644444E ♦O 1 
.62)6687E^01 
,632POOOE^01 
,640P337E^0l 
.6500000E>01 
,P87F471E^01 


.?09"333E^O? 
.?02-’333P*0? 
.?22R330E^0? 
.?22qSP0F^0? 
,?09n33lE^02 
.?R614?9F^02 
.3047333E»02 
.3E53000E^02 
,378'’333E*02 
,3877273'^*02 
,389ROPOE^02 
,3P6(t6f  7F»02 
,385(^Q23E^02 
.374''OOOF^02 
.366P1P2P»02 
,3??pQ47F*02 
.72n905F^0? 
.3P9p?2?E^02 
.2P3'==000E*02 
.?75«000E*0? 
.?52->P46F*02 
,239r'OOOE*0? 
,167ROOOE»0? 
,l64^ft00F*02 
.1680POOF4.02 
.1711667F402 
,)71?IP5F*02 
.177PPP7E*P? 
.1801)  U'^^02 
.1P21P79F*0? 
,190aOOOE^O2 
.22?f^8P0E^02 
,21?q4l2E^02 
,?07t333E»02 
.206???2E»02 
,?05«-667E^02 
,209'^P33E^02 
.?12'^0ft0E*02 
.213'^'Q09F*0? 
.21?f’4l2E*P? 


L. 
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C0^  TOUR 


.7000^*01 
.7000F+01 
.7!)00F-f01 
.7000F*01 
,7ono'^*oi 
.7000F*01 
.7000F+01 
.7000F4.01 
.7000F+01 
,7000F*01 
.7000F+01 
.7000F*01 
.700oe^oi 
.7000F*01 
,70ft0F>01 
.7000‘^  + 01 
,70OOF*01 
.7000F+01 
,P000F*0l 
.flOOOF*01 
,fi000E*01 
.R000»^^01 
.POOOF+01 
.ROOOf *01 
.P000F4.01 


Y 


,?^PR0RTF^0? 
.24PO?77F>n? 
,?1F437FF^07 
,?1?OOOOF*0? 
,2?1705RE*0? 
,?71 12S0E^OR 
,?6^5000E*OP 
.306P1P?E*0? 
,33?4737E*0? 
,357PF71E^0? 
,3R7SOOOF^O? 
.366Plfl?F^02 
.33=^?S00E*t)? 
,3070R0PE^02 
,300F667E>0? 
,?4<?lP4qE^n? 

.24PQOR3E^O? 

.?SF3(S3PF^0? 

.25?3l2SE^ft? 

,?5?0000F*0? 

.?5?««24F>0? 

,?5737S0E*0? 

,?660000F^O? 

.255363FF4.0R 


.?6?PR30F*0? 
.?63'^2P4F*0? 
,?P1 1R75F+0? 
.343637SF+0? 
,3471  17E.E1-0? 
,345‘=:000F^0? 
,33P='333E^0? 
,3304S45E*02 
.3127F.84E^02 
,3021429F»02 
,2617C00E*02 
.231«1P2E^02 
,?197500E^02 
,226P091E»0? 
.234?22?E^02 
.262o291F*02 
,2P2P720F*02 
.?62Q930E^02 
,303^000F^02 
,305'=:6?5F*0? 
.31071?5F*02 
,310'^471F»02 
.31C<^000F*02 
,306'^OOOF  + 02 
,303r.0O0F40? 


Appendix  D. 
PROGRAM  LISTING 
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noo  n n too 


OVFPLAY (CONTOUR, O.n) 

PROGRAM  MAIN(TAdf1=65,0IJTDUT=6'"  , INPUT  = TAPF1  , T aPF5  = ^S  , T A0F6  = 65  , 

1 TAPF7  = sn,TAPF?=S13,TAPF9Q) 

THIS  CONTOIIO  PROORAN  DFVFLOPEr)  “Y  OPAOY  PATRICK  — FFR  1R77  — 

AT  MIRAOCOM  REOSTONF  ARSENAL 

OIMFNSION  NL(T) 

C0MM0N/C0NTRL/NC0NTR,VAL (?0) , KLOW, INM, JMAX.NFLM, IPT 
C0MM0N/mxmn/XM4X,XMIN. YMAX,YMU’, XC, YC 
CONMON/HIdloT/ImI , ITK 
cpmmon/corf/kopf 
MCONTP=0 
IFLOW=l 
INM=?0 

JMAX=n 
NFLM-n 
IPT  = A 
IHI=0 
ITK  = 0 
KnRP=0 
REWIND  I 
REWIND  ?. 

REWIND  7 

THESE  STATEMENTS  ARE  FOR  SYSTE^'  LOADINS  IN  THF  MAIN  OVERLAY  ? 

TF(.|MAX  .1^0,  (!)  ort  TO  10 
RACKSPACE  7 
PEAD(7)  NFlM 
10  CALI.  COMNFC(S) 

CALI  COMNFC(6) 

WRITE (ft, 9700) 

FSTARLISH  roNTROL  PARAMETERS 
70  CALL  TEK(1) 

WRITE  (6,91 00)  TFL0w,IMM,IPT,IHT,ITK,NC0NTR 
READ(S,R000)  NAme,IVAL 
IE(NAmF  .EO,  AHIFLO)  IFL0W=IVA| 

IE(NAME  ,Fo.  AHINiim)  INM=IVAL 
IF(NAme  .FO,  AhIPTS)  IPT=IVAL 
IF(NAME  ,fo,  ahIhtp)  ihi=I 
IE (NAME  .FO.  AMTTFK)  DO  TO  30 
IF (NAME  ,FO.  AHNCON)  NC0NTR=IV»L 
IE (NAME  ,F0.  4HCONT)  on  TO  AC 
IF (NAME  ,FO.  3HFNO)  GO  TO  SOO 
GO  TO  ?0 

INITIALI7E  TEKTRONIX 

yO  ITK=1 

CALL  INITT(O) 

CALL  TEPM(3,107A) 

CALL  CHRSI7(3) 

GO  TO  ?0 
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40  IFdHI  .NF.  1)  r,0  TO  SO 
REWIND  RR 

CALL  INITAL  (RRtno.2?.0) 

CALL  RSTRd) 

SO  CONTIMUF 

IFdFLOW  ,FO.  1)  r.O  TO  300 
CALL  TEKd) 

WRITE (6«R?nO)  NCOMTR 
REAO(S»»)  (VAl (N) .N=l .NCONTR) 

IFdFLOW  ,RF.  P)  oo  TO  300 
RFAOd,*)  JMAX.MEL'^ 

WRITE(7)  IMAX.NFLm 
DO  100  I=I , JMAX 
REAPd,*)  X,Y»7 
100  WPITF(7)  X,Y«7 

IF(MELM  .LF,  0)  on  TO  300 

DO  ?00  1=1 .NELM 

PEADdt*)  NL  d ) »NL  (?)  .NL  (3) 

?n0  WRITE(7)  (NL (K) ,K=) ,3) 

300  JFLOW=IAas(IFLOW) 

DO  400  I=JFLOW,3 

CALL  OVFRLAY(7HCONTOUPd  <0) 

400  CONTINUE 

CALL  TEK(l) 

IFLOW=3 

IVAL1=IVAL 

IFdHI  .GT.  0)  CALL  RSTR(?) 

GO  TO  10 
SOO  CONTINUE 

IFdHI  .GT.  0)  CALL  RSTR(?) 

WRITE (6.R400) 

CALL  DISC0M(6) 

CALL  EXIT 
flOOO  F0P^'AT(A4,I?) 

RlOO  FORHATdH  * NAME*7X*VALUE * 1 OX»FORMAT  A4d?*/ 


A 

R 

C 

0 

F 

F 

G 

H 

I 

J 


* IFL0*RXI2,4X*HINUS  to  OMIT  POINT  NUMBERS*/ 
1RX*1*4X*0NLY  X.Y  ANp  7 COORPINATFS  APE  INPUT*/ 
;SX*?*4X*C00PDINATFS  AND  ELEMENTS  ARE  iNRUT*/ 
1SX*3*4X*RERL0TTING  PF  GFNERATEP  CONTOUPS*/ 

ff  INt)M*oxI2.4X*SPACI^'G  OF  CONTOUR  NUMBERING*/ 

* IPTS*10XI1 »4X*SM00THING  FUNCTION  FOR  SPLINF  FIT*/ 

* IHIP*10XI1 »4X*SFT  TO  I FOR  HOUSTON  INSTRUMENT  PLOTS*/ 

* ITFK*10XI1 »4X*SET  TO  1 FOR  TEKTRONIX  OUTPUT*/ 

* NCON*10XI1«4X*NUMPER  of  contours  to  be  plotted  20MAX*/ 

* CONT*1 0X*0*4X*TO  CONTINUE*/ 


K * END  *I0X*0*4X*TO  E^'D*//) 

9200  FOPMATdH  * ENTER*I5*  CONTOURS  IN  ASSEnDING  ORDER  FRFF  FIFl  D*/) 
R300  FORNATdHO) 

R400  FORMAT(IHR) 

END 
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lOENT 

kOREEL 

ENTRY 

KORFEL 

EXT 

CPC 

LIRT 

D 

TALI 

KORfFL ( 

A,  K(1RE.  KFI  ) 

KFL  = 

FIELD  length 

KORE  = 

XEL  - ADGRESG  of  A-1. 

VEO 

•»f)/nHKDREEL  «?A/? 

AO^AVr 

HSS7 

1 

KOOTFL 

nATA 

0 

SX«S 

AO 

<;pl 

XI 

FTN  CONVERSION 

«:A1 

A1  *1 

ETN  CONV'^P'^ION 

FA<> 

AORAWE 

«:R? 

XI 

ETN  CONVERSION 

RAl 

Al  ♦! 

ETN  CONVERSION 

RRT 

XI 

ETN  conversion 

SXib 

R1 

FA*. 

HI  RAVE 

RAVE  Rl  (ADDRESS  OF  A) 

qxF 

p? 

RA6 

R?RAVE 

SAVE  H?  (AODPESS  OF  KORE) 

RX^ 

HI 

SAfi 

R3RAVE 

RAVE  H3  (ADDRESS  OF  KEL) 

MX6 

0 

RA^ 

length 

Memory 

CM, length, RECALL 

7EM0RY  MACRO  CALL 

RAl 

length 

LOAD  ( ENGTH  INTO  xl  (((PPro  HALE) 

l.xi 

30 

RHIEJ  INTO  LOWER  male 

RXft 

Xl 

field  LENGHT  into  X6 

RA? 

B3RAVE 

SAft 

%2 

STORE  LFMGHT  into  USER  WORD 

RA7 

RIRAVE 

addrfer  of  a into  X? 

TXR 

Xl-X? 

RX^ 

X3-1 

KORE  INTO  X6 

RAl 

P?RAVE 

RA6 

XI 

STORE  KORE  INTO  USER  WORD 

RAl 

AORAVE 

RAO 

XI 

EO 

KORFEL 

LEMATH 

DATA 

0 

storage  for  LENGTH  PETtlRNFO  HY  C 

R1  «;avf 

nATA 

0 

storage  for  hi 

n?cAVE 

nATA 

0 

storage  for  b? 

R3«^AVE 

nATA 

FNn 

0 

STORAGE  EOR  HI 

'1 


OVFPLAY (CONTOUP.l .0) 
oonr,RAM  TPIAN^' 

rnMMON/rOMTRL/NrOKiTR, VAl  (?0)  , INUM.  fMAXtMFl.M,  IPTP 

ro«MON/rORF/KOPP 

common  A(U 

CALI  KOPE^L (A,KOor,KrL ) 

KOP'^  = ‘iA*  (KnRF/E4) 

no  IOC  T=»,KOPc 

100  A ( 1 ) =0.0 

PEWTNO  7 

PEAn<7)  JMAX.NFI” 

J4=0*JMAX 
LK-5  IMAX*  1 
L^l =2*J“AX*1 
Lx?=3*  IMAX  ♦ 1 
Lk7=4»JMAx+1 
L>4  = S*  J^^AX  ♦ 1 
Li<''i=6*./^'AX«  I 
MK= AOJA 
MK  = MK  ♦l.K=; 

LOJ(-  = k0PE-mk 

li^lLOir  .TT.  0)  GO  TO  ?00 
CALL  TFK(!> 

WPTTF (ft.POOO)  mk.kOPF 
CAI  I EXIT 

?O0  fail  TPI(A(1),A(L'<),A(lk1).A(L''?).A(LK1),A(LK4),A(IK5).JMAX.NFL*') 
QOOO  ropMATdHI*  COPF  TNSUFFICIFNT  I FEDFO  =«17o  AVARARIF  =»I7) 

FNO 
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<;UPO0UTINr  TRI  ( X . v , 7 , xP . YR . I OM I T . I FL  M , jMA x , N 1 ) 


1 


ThIF  ROI.ITI^F  PORMIILATPS  TRTANjGi  FS  FROM  JHF  NODAL  POINT  DATS 

ni^'FNSION  X ( 1 ) «Y(  1 ) *7(1  ) «XR  ( 1 ) , YR  ( 1 ) , TOMIT  ( 1 > . IrLMCI,  1 ) , n (S) 

?0  DO  110  IOrl,JMAX 
ilO  IOMIT(IO)=0 

CALI  RFAD(X,Y,7, JMAX) 

OFTFRMINE  THE  CENTER  OF  THE  NODAL  POINT  DATA  TO  BFOIN  THE  TRlANOLFS 

CALL  CENTER  (X  .Y.J-^AX.  Jl  ) 

1?0  CALL  CL0Sr(X.Y*7.  twAX,Jl,J?,j3.IOMlT,l) 

CALI.  ROT(X,Y.XR.YR,  JMaX.  J1 
CALI  CLOSF  (XP.YR.7,  JMAX.  J1  *.)?*  )3.T0MTT.?I 
TP  ( 13  ,NE.  0)  00  TO  UO 
lOMIT  ( J1  )=:0 
J?  = >2 
00  TO  1?0 
160  CONTINUP 
N=1 
Nl  = l 

TFLM(l,Nll=Jl 
IFL«(’.NH=J 
IPl  M n.Nl  ) = J3 
100  IL(1 »=IPLM(1.N1 ) 

IL (?)=IPLM(?.N1 ) 

11  (3)=lPLMn.Ml) 

IL  (6)  =1FL‘'  ( 1 *N)  ) 

IL (P) =IFLm (?,N1 ) 
no  500  M=i ,3 

CALL  compare (N.Nl . IL (w) , II (M. 1 1 , IPLM, irOM) 

IP  (I  COM  ,r,T.  0)  00  TO  POO 

CALL  POT (X.Y.XP.YP.JMAX.IL(M).Tl<M.l).TL(M.P)) 

n=o 

IfO  CALI  CLOSP (XR.YR.7, JMAX.IL (M) .TL (M.n  . I3.I0MIT.?) 

IP(  13  ,P0.  0)  on  TO  ?0O 

CAI  I oinFOK  (N.  IL  (M)  . IL  (M.n  . J3.  IFLM.  IChK) 

IPdCHK  .PO.  0)  00  TO  1(SS 
I0MIT( J3)=? 
on  TO  lAO 

CAI  L ChPCk ( J3. IPLM.X. Y.N. ICHk . tOMI T) 

IPdCHK  .=-0.  0)  00  TO  170 
IDMIT ( J3) =1 
00  TO  160 

170  IP(J3  .FO.  0)  00  TO  ?00 
N = N.  1 

IFLMd  .N)  = IL  (M) 

IPLM (?,N) rTL (M.l ) 

IELM(3.N)=.n 
700  CONTINUP 

IP(N  .FO.  Nl)  00  TO  300 
N1=N1*1 
00  TO  100 
300  COMTTMIIF- 
REWIND  7 
WPTTE(7)  JMAX.Nl 
no  600  J=1 .JMAX 
600  WPTTP(7)  » ( J) .Y( J) .7( J) 

DO  SOO  N=1 .MI 

SOO  WPITF(7)  dPLM(K.M) .K=l .3) 

POD 


J 


r>  o o D 


SUPBOUTINF  REAn(X.y.Z, JMAX» 


THIS  ROUTINE  READS  THr  INPllT  DATA  AND  FINDS  The  MAX  AND  MtM 
AMD  THE  COORO  OF  THE  TENTER  OF  THE  PLOT 

COA'MON/MXMN/XMAX.XMIN,YMAX,YMl*i,XC»  YC 
dimension  X(1).Y(1).Z(1) 

XMAX=-1 ,E?0 
YMAX=-I .E?P 
XMIN=  l.E?0 
YMIn=  1,E?0 
IFi  AG=0 

DD  ion  J=1,JMAX 
RFAD(7)  X ( J) fY ( J» ,7 ( J) 

ino  lF(x(J)  .t  T.  n.p  .OR.  Y(j)  ,lt.  o.O)  iflas=i 
IF(IFLAG  .FO,  1)  T.O  TO  500 
ZOO  DO  TOO  J=1 , JMAX 

IF(X(J>  ,r,T.  XMAX)  XMAX  = X(J) 

IF(X(J)  ,LT.  XMIN)  XMIM=X(J) 

1F(Y(J)  .GT.  YMAX)  YMAX  = Y(,J) 

1F(Y(J)  .IT.  YmTM)  YMTN=Y(J) 
mo  rONTINUF 

XT= (XmaX*XmIN)*O.F 
YT= ( YMAX*YMIN) *0.S 
RFTIIRM 
qOO  XM=0.0 
YM=0.0 

DO  GOO  .1=1  . JMAX 
1F(X(J)  .LT.  XM)  XMrX(J) 

IF(V(J)  .1  T.  YM)  YM=Y(J) 

GOO  CONTINIJF 

DD  700  J=1 .JMAX 
X ( J) =X ( Jl -XM 
Y ( J) =Y ( J) -YM 
700  CONTINUE 
GO  TO  ZOO 
END 
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onn  OOO  00*^ 


SUROOUTINF  CLnSF(<.Y»7tJMAX.Jl . J?tJ3f InMITtKl ) 


C 

C THIS  ROUTIK'F  FINOS  BASFO  ON  Kl 

r Kl  = 1 THF  ? POINTS  NFAPFST  TO  jj 

0 Kl  = 2 TWF  DOINT  NFAPFST  THF  AVEpAGF  OF  J1  AN'O  j? 

r 

niMCNSioN  X (1 ) ,Y(i » «zn ). lOHiTn ) 
COHMON/^*XHN/XHAX,XHIN,YMAX.YMIm,XC»YC 
“=0  SCALF  = 0.0 

XLNIT=1.2*AHS(X(J1 ) ) 

J3=0 

XP=0.0 

YP=0.0 

IFIKI  .MF.  1)  on  TO  100 
J?=0 

XP=X ( J1 ) 

YP=Y ( J1 > 

100  R1=1.F20 

200  00  SOO  J=1 , JMAX 

IFdOHITlJ)  .GT.  0)  GO  TO  SOO 

1F(J  ,F0,  J1  ,0P.  J ,F0.  JP  .Od.  J ,F0.  J3)  GO  TO  SOO 
IF (Kl  .FO.  1 ) GO  TO  PRC 

CHFFK  TO  SFF  IF  POINT  IS  WITHIm  XLIMIT 

IF(ABS(X(  in  .LT.  XLMTT  .ANO.  Y(,()  .FtT.  SCALE)  GO  TO  290 

CHPCK  IF  POINT  IS  GREATER  THAN  SCALE 

IF(V(J)  .LT.  SCALE)  GO  TO  SOO 
JJ=  )1 

CHFCK  IF  POINT  IS  WITHIN  45  OE G OF  J1  ANO  J2  LINE 
A = 1 .0 


IF  (X  ( J) 

.1  T. 

('.0 

.ANO. 

X ( J1  ) 

.1  T. 

0.0) 

A=-l .0 

IF (X ( J) 

.I,T. 

0.0 

.AND. 

X ( J?) 

.1  T. 

0.0) 

GO  TO  220 

(X  ( J) 

.gt. 

0.0 

.AND. 

X ( J1  ) 

.LT. 

0.0) 

JJ=J2 

2)0  XX  = A*(X(J  ))-X(J) ) 

IF (XX  .GF.  0.0)  GO  TO  290 
IF (Y ( J) /ADS(XX)  .LT.  1.0)  GO  TO  SOO 
GO  TO  290 
PPO  JJ=,)2 
A=-l .0 
GO  TO  210 

5Q0  R=SORT  ((X(.I)-XP)**2*  ( Y ( J)  - YP)  **2) 
IF(P  .GT.  Pi)  GO  TO  SOO 
PJ  =P 

GO  TO (300.400) «K1 
300  J?=J 

GO  TO  500 
400  J3=J 
500  CONTINUF 
500  CONTINUF 
RETURN 
ENO 
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SlJCiPOUTTNF  ROT(x,Y»XR,YR.  JmAX.  II . J?.  Jl) 


***  THI';  RODTINF  ROTATFS  ALL  POINTS  SO  THAT  Thf  X AXES  LIFS  A|  ONG 
***  JI  TO  J?  AMD  THF  Y AXFS  IS  OIPrCTFO  AWAY  FCO^'  JO 
»**  WITH  X ( Jl ) = -X ( I?) 

DIMFNSION  X (1 ) ,Y (1 ) « VP(1 ) ,xP( 1 ) 

0=X ( J?) -X  ( Jl ) 

IF(D  ,NE.  0.0)  GO  TO  ?00 
SM=1 .0 
CS=0.0 
xiNr=x (Jl ) 

GO  TO  300 

30  0 IF(Y(J3)  .MF.  Y(JD)  GO  TO  ??0 
X‘^=0.0 

DO  ?10  J=1 . JMfiX 
XP( J)=X(J) 

YR( J)=Y( J)-Y( Jl ) 

310  GONTIMUF 
GO  TO  110 

??0  Xw=(Y(J3)-Y(J1 ) ) /n 
THFTA  rATAM(XM) 

SN=FIN(THFTA) 

CS=COS(THFTA) 

3F0  XINr  = X ( Jl ) -Y ( Jl ) /XM 
300  CONTINIIF 

DO  100  J=1 . JMAX 
XP  t J) =X ( J) *rS  + Y ( J) «SN 
100  YP ( J) =XING*SN*Y ( J) *CS-X ( J) oSN 

SHIFT  X AXFS  TO  THE  CFMTFR  OF  ll  AMD  Jp 

iio  XAVG=  ( XP  ( )1  ) *XP  ( 1^)  ) *0.5 
DO  ISO  J=1 .JMAX 
ICO  XP  ( I)  =XP  ( J)  -XAVG 

IF(J3  .FO.  0)  GO  TO  SnO 
IF(YP(J3)  ,lE.  0.0)  GO  TO  coo 
DO  AOO  J=1.JMAX 
AOO  YP  ( I)  =-YR  ( J) 

SOO  COMTINDF 
RFTDRM 
Fun 


I 
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SUPOOUTINF  CHFCK ( JT . I , x . Y ,N : , TCHK, IomIt) 


1! 

r 


***  THIS  OOUTIMF  CHFCKS  To  SFF  TF  THF  JOINT  THAT  H4S  RFFN  FOUMO 
DOFS  NOT  HAVE  ALL  THE  pOSS I « I L I T I FS  USFD 

n I VPNS  I ON  X(l),Ym»lFLM(3.n«T0HITm 
IF(I0VIT(,)T)  ,r,T.  0)  00  TO  ?0OO 
PI=3.1415Q?^ 

ICHK=o 
ANGl=0.n 
ANO=P  .0 

DO  1000  K=1 ,N) 
no  100  1=1.3 

IFdELMd.x-)  .FO.  JT)  00  TO  ?0'' 

IPO  CONTINUE 
on  TO  IPO" 

POP  on  TO  (?io.??n.?3r>)  .1 

21P  J1=IFLM(2,k) 

J2=IFLM(3,K) 
on  TO  300 
?PP  J1=TFLM(1,k) 
jp=IELM(3.k) 

00  YO  300 
’30  Jl  = IEL'^d.K) 
jr>=IFLV(2,K) 

300  A’= (X ( J1 ) -X ( J2)l «»?♦ ( Y ( Jl ) -Y( J2) ) 

P = SORT  ( (X  (,)T»  - X (J1  ) ) »«?♦  ( Y ( )T)  - Y ( J1  ) ) ««2) 

C=SORT ( (X ( JT) -X ( JP) ) **2* (Y ( JT) -Y (J2) ) “02) 

0= ( A2-R»B-C*C) / (-’.*R*C) 

IF(ARS(0)  .LE.  0.Q909PR)  00  TO  400 
AMO) =0 .0 

IS(n  ,LT.  P.O)  ANr,l=Pl 
on  TO  ROO 
400  rOF'TINUE 
SION-1  . 

lF(n  ,LT.  0.0)  SION=-l. 

A^|01=ACnS(ARS(|■))  ) 

If (SION  .OT.  0.0)  00  TO  ROO 
ANM=P(-AM01 
ROO  CONTINUE 

AM0=AN0+ AMOl 
lOnO  CONTINUE 

ThET=2.*PI-AN0 
IF(THST  .LS.  n.POl)  I0HK=1 
’OOP  RETURN 

END  ' 
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SUgPOUTINr  COMPARF (N.Nl . J1 ♦ J?, TFLM, ICOM) 


»»*  THIS  ROUTINE  CHECKS  TO  SEE  IF  THF  JOINTS  J1  AND  J2  ARE  FOUND 
***  IN  ANY  OTHER  TRIANGLE 

dimension  IELM(3»1) 

IC0M=0 

DO  400  K=1 ,N 

IF(K  .FO.  Nl)  GO  TO  400 

no  100  Kl=l,3 

IF(J1  ,NE.  IEI.M(K1,K))  GO  TO  I’^O 
GO  TO  ?00 
100  CONTINUE 
GO  TO  400 
300  no  300  Kl  = l ,3 

IF(J?  .NF.  TELM(K1,K))  GO  TO  3''0 
GO  TO  500 
300  rONTINUE 
400  CONTINUE 
RETURN 
500  1C0"=1 
PE  TURN 
ENO 
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I 
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f 

I 
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stiBPOUTTNr  sinrrK  (n,  Ji , J2» n, i^lm.  ichk i 

c 

r THIS  suPRoiiTiNir  rnrcKs  to  sfs  tf  a psoonsFo  tpianoif  pops  mot 
c HAVF  A SIPF  THF  same  AS  TWO  OT^iFP  TPIAmOLFs 

r IT  OFTUPNC 

r TTHK  r 0 NONE  FOi'NlO 

r TTHK  = 1 T'VO  sro^s  FOONO 

r 

n I MENS  I ON  tflm  (T,  I ) ,j  (0) , jcHK  ('!).iin).i2n) 

DATA  12/2, T»l/ 

DATA  11/1,2,3/ 

00  ion  JJ=1,3 
inn  jCHK(jj)=o 

J(l)=Jl 
J(?) =J2 
J(3)=J3 

c 

C LOOP  ON  3 POSSIRII  ITIES  IN  A IpTANOLF 

r 

DO  FOO  Kl  = 1 ,3 
Ml=l 1 (K1 ) 

M2=12  (Kl ) 

LOO'’  ON  ALL  elements 
no  SOP  K=1 ,N 

Chfck  joint  J(mi)  on  all  joints  IN  ELEMENT (K) 
no  200  M = i ,1 

IF()(M1)  ,FO,  TFLM(M,k-))  GO  TO  onn 

POO  rOMTIMUE 
GO  TO  SOO 

IS  ONE  IS  FOUND  OHFCk'  IOINT  J(M2)  AN  ALL  JOINTS  IN  ELFmEnt(K) 
TOO  no  aoo  m=i  ,1 

1F(J(M2)  .FO.  TFLM(m,k))  JOHk  (k-J)  = JCHK  (Kl)  *1 
AOn  COMTINUF 
SCO  CONTINUF 
FOO  CONTINUF 
C 

r CHFCK  IF  TWO  STOPS  WEPF  FOUND 

C 

1 GMK-xn 

DO  ■'00  JJ=1  ,3 
IF  ( ICHK  ( J ))  .GT.  1 ) ICHK  = 1 
700  CONTINUF 
RF  TURN 
FNP 
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OVFPLAY (COMT0UR.?.0) 

PPOr.RAM  mpcOMTR 

COMMON/rON'TRL/r.'COMTR.  VAL  (?0)  , I FLOW , I MJm  , JMfl  x , NELM , IPTS 

COMMON/rOPF/KORF 

COMMON  A ( n 

CALL  KuPEFL (AtKORF.KFL ) 

K()RF  = 64*  (KORE/64) 

00  100  T=1.K0PF 
ino  A (I)=0.0 
PFWTNO  7 

RFAn(7)  JMflX,ME|.M 

LK=JMAX*1 

LKl=2*JMAx*l 

LK5=3*jmax*1 

mk=o*NFlm+lk? 

LOIF=KOPE-mk 

IF(L0IF  ,0T.  0)  GO  TO  AOO 
call  C0NNFC(6) 

CALL  TEKM) 

WRITE (FuROOO)  LOTF.KORF 
CALL  EXIT 

AOO  CALL  CONT (A ( n ,A  (LX)  • A (LKl ) , A C.K?) ) 

no  GOO  1=1 ,KOPF 
GOO  A(l)=0.0 

-kg=kope/s 

LK1=2*KG*1 
LKP=4*KG+ 1 
CALL  0IFCnM(6) 

CAl  L ORDER  (A  (n.A(LKl),A(LK2)  .k-S.  INUM) 

C«IL  C0NNFr(6) 

ROOD  FOPMATUH  * CORE  INSUFFICIENT  i''EFnFO  =*17*  AVAILABI  E =*I7) 
E^in 


.i 


t 


iJ 
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SURROUTINF  COWT ( X , Y f Z . IFLM ) 

DIMFNSION  X(1)  tY(l)  «Z(1) 

COMMON/CONTRL/NCONTR.VAL  (?f>)  , I FLOW  , I NU^  , JM  AX  .NFLM  , T PTS 

A CONTODRTMG  PROGRAM  rOR  RANDO'''  SPACER  POINTS 
THE  POINT*;  MAY  RF  RANOOMLY  NUMRFPEO 

ALL  POINTS  MOST  RF  CONNECTED  TN  SUCH  A MANNER  TO  FORM 
NON-OVEPLAOING  TRIANGLES  CALLFO  ELEMENTS 

CALL  TEK(l) 

READ  THE  X - Y coordinates  AND  THE  NOoAL  VALUE 

POO  DO  POO  J=1 . JMAX 

POO  RFAO(7)  X(J) ,Y(J) »7(J) 

READ  THE  NODES  FOR  EACH  ELEMENT 

00  AOO  N=1 ,NELM 
400  REAO<7)  < IFLM(K,N) «K=1 ,3) 

REWIND  7 

WRTTE(7)  IMAX.NFLM 
DO  410  J=I.JMAX 
410  WRITF<7»  * f J) f Y f J) .7 ( J) 

DO  420  N=1 ,NELM 
4P0  WRITE(7)  ( IELM (K ,N) «K=1 »3) 

CALL  MAXMTN(XtY,JMAX.TFLOW) 

CALL  C0NT0UR(XfY.7,IELM) 

SOO  CALI  DISCON(6) 

WRITE (6»9000)  NCONTRt JMAX,NELM 
WPITE(6«9200) 

DO  ISOO  1 = 1 t JMAX 

600  WRITE (6*9400)  I * X ( I ) , Y ( I ) . 7 ( 1 ) 

WRITE (6.9700) 

DO  700  L=1.NELM 

700  WRITE(6.9S00)  L . IFLM ( ] ,L ) ♦ I ELM ( 2 *L ) . lELM ( P *L ) 

REWIND  2 

ROO  REA0(2)  K,A,H,I 

IF(F0F(PM  1300,900 
ROO  IF(K  ,GT.  ?)  GO  TO  1200 
GO  TO  (1000,1100)  ,K 
1000  WRITE (6.9P0O)  T.A.R 
GO  TO  ROO 

1100  write (6,9000)  A,R 
GO  TO  800 

1200  WRITE (6,9600)  I, A 
WRITE(6,9700) 

GO  TO  800 
IPOO  CONTINUE 

CALL  C0NNFC(5) 

CALL  CONNFC(6) 
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9000  FOPMATdHl  *Niimpfr  OF  CONTOURS  TO  SF  PlOTTFO  =*15/ 

1 * NIJMPFR  OF  NODES  =*15/ 

? * niimREP  of  elements  =*I5/) 

9100  FORMATdHO*  CONTOURS  TO  PE  PLOTTED  *//) 

9200  FOR’^ATdHO*  NODE*A  X*X-COORD*l  1 X*  Y-C00RD*1TX*V  ALUE*// ) 

9300  FOOMATdHI*  element  NODE  1 NODE  ? NODE  3*//) 

9400  F0RMAT<I5,^E1P.7) 

9S00  F0RMAT(4IO) 

9600  FORMAT dH|*CONTOUP  NUMBER  *T5*  CONTOUR  VALUE  =*E14.7) 

9700  FORMATdHO*  ELEmFNT*29X*L 1 NE*/1 2X*X* 14X*Y*1 9X*X* 14X*Y*//) 
9S00  F0PMAT(I5,?EIP.P) 

9900  FORMAT (1H*,39X,2E1S. 81 
END 
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SURROUTIN'^  COMTOI)R(X«Y»Z«IFLM) 

COM»<ONI/COMTRL/MCOMTR«  VAL  (?0  ) » 1 FLOW  ♦ I NIJm  , ju  a x *NFl  M,  TPTF 
CO'-'”ON/MIplOT/IHI  ,ITEK 

niMFNSION  X(l),Yn)fIF|Mn,l),r(4>,TK(A).TY(<.) 
niMFNSION  LABL (?0) . IAR(6) ,7( 1 ) 
no  100  1=1 .?o 
100  LARL(I)=0 

DO  1100  K=1 .MCOMTP 
N = 0 

NCONTRl=NrnNTR+10 

WRITE  (?)  MCONTRl ,VAL (K) .VAL  (K) .K 

no  1000  J=1*NFIM 

J1=TELM(1,J) 

J?=TELM(2, J) 

JO=IELM(3, J) 

F(1)=Z(J1) 

F (?) =7 ( J2) 

F (0)=Z( Jl) 

F(4)=F(1) 

TX(1)=X(J1) 

TX(R)=X(JR) 

TX(0)=X(JO) 

TX (4)=TX (1 ) 

TY ( 1 ) =Y ( J1 ) 

TY (?) =Y ( J?) 

TYn)=Y(J?) 

TY(4)=TY(1 ) 

IF(K  .r,T.  1)  r,0  TO  100 
IFdTFK  .MF.  1)  no  TO  10  0 
CAI  L MOVEA  (TX  ( 1 ) «TY  ( 1 ) ) 

no  ?oo  M=?,4 

?00  CALL  ORAWA (TX (M) ,tY (M) ) 

100  CONTINUF 

BA  = AMIN1  C^(l)  .F(?)  ,F(1) ) 

AA  = AMAX1  (r (1 ) ,F (?)  ,F (1) ) 

IF(VAL(1)  .OT.  AA  .OR.  VAL (MCO-ITR)  .LT.  BA)  CO  TO  IlPO 

JF(VAL(K')  .CT.  AA  .OR.  VAL(K)  .LT.  BA)  CO  TO  1000 

IRCM=0 

XX1=0.0 

YY1=0.0 

no  ROO  L = 1 .3 

T''(F(L)  .RO.  VALC<))  F(L)=F(L)*1  .0001 
1F(F(L*1)  .FQ.  V/AL(K))  F(L*1)=F(L*1  )*1.0001 
IF  ( (F  (L)-VAL  (K)  ) * (V/AL  (K)-F  (L*l  ) ) ) ROO.BOO.BOO 
con  ALP=  (VAL  0<)  -F  (I.)  ) /(F  (L*!  ) -F  (L  ) ) 

XX=TX (L) *ALR* (TX (L*l ) -TX (L) ) 

YY=TY (L) ♦ALP*(TY(L*1 )-TY(L) ) 

IFdRCNI  .MF.  n r,0  TO  ROO 

IF(xxi  .Eo.  XX  .A'jn.  yyi  .eo.  yy)  co  to  boo 

IFdTEK  .FO.  1)  CALL  OABHA  (XX  .YY.1) 

IRCM1=IRGM*1 
WRITE(?)  IRON! ,XX,YY,J 
LARL  (K) =LARL  (K) ♦! 

I B=LARL  (K) 

IF  (MOn  (LR.  INU^*)  ) 700«RO0»7O0 
COO  IFdTFK  .Mr.  1)  CO  TO  700 
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CA| L CMDS17(4) 

CALL  ffO»H(vAL  (K  ) .(t.,  1 , 1 ftP.l?) 

CAU  JUSTPD(6, I AB.O»3?.LFN*TOF=ET) 
TALL  ‘'OVPFL  ( lortiET.O) 

CALL  HLARPL  (I  ) ) 

700  IQr,*c  = 0 
N=r 

r,o  TO  POO 

OOO  IFCITTK  .FO.  n CALL  wnVEA(XX.YV) 

IF(^I  ,NF.  0)  PACKSPACF  ? 

WRTTF(?)  inr,Nl  ,XX,YY.  J 
XXI =XX 
YYl =YY 

IHRN=1 

M=1 

900  COMTImljf 
1000  inc*l=0 
1100  CONTINUF 

IFdTEK  .ME.  1)  RFTUPM 
CALL  BELL 
CALL  TST'Nn 
CALL  TTMPMT(ICP) 

RETURN 

EMH 
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SURROUTINF  MAXMJ N ( X ♦ Y , JMAX , TFL ) 
COMMON/MXMN/XMAX.XMIN,YMAXf YMlN.XCt YC 
COmmON/HIPLOT/ImI , itfk 
niMFNsinM  xn),Y(i) 

IFdAPSdPLOW)  .FO.  1)  GO  TO  ?t0 
XMAX=-1 .EPn 
XMTN=  l.E?0 
YMAX  = -1  .EPr> 

YMTN=  l.E?0 
no  100  J=1 . JMAX 
XMAX=AMAX1 (XMAX.X ( J) ) 

XMIM=AMIN1  (XHIN.X  (J) ) 

Ymax=AMAX1 ( YMAX«Y ( J) ) 

YM1N=AMTN1 ( YMIN.Y ( J) ) 

100  COMTIMUF 

?00  IFdTEK  .ME,  1)  GO  TO  TOO 
CAIL  SWINnn(200,6S0.S0.650) 

CA|  1.  OWINnn(XMIN,XMAX,YMlN.YMAX) 

CALL  TFK(l) 

300  CONTINUE 
RETURN 
ENO 
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SURROUTTNP  OROFR ( X . Y * TFL AG .LMA * . I NUM ) 

THIS  ROHTI^F  ORDERS  THE  CONTOUd  LINES  TO  EOPH  A CONTINOUS  I INE 

CPMMON/HIPLOT/IHI , ITEK 
DIMENSION  X(2.n  fV(2fl) ,IELAG(1  ) 

REWIND  ? 

ISTOP=0 

inO  REAO(2)  K .CONtCONl  ,KK 
IE(EOE(?))  2<»00,?00 
200  MORE=0 

IE(t<  .GT.  100)  GO  TO  2A00 
DO  TOO  I=1,LMAX 
300  IFLAG(I)=-I 
CONT=CON 
N=1 

IPEO=0 

400  READ(2)  Jl.A»o«KK 
TE(POE(?))  2S00.500 
500  CONTINUE 

IE(  ll  .GT.  ?)  GO  TO  900 
X ( J1 ,N)=A 

Y(J1.M)=B  1 

READI2)  J?.C«n,KK 

IE(EOE(?)>  2400.600  i 

600  CONTINUE  ■ 

IE(J?.GT.  2)  GO  TO  900  ! 

x(j2.N)=c  ; 

Y(J2.N)=0  S 

71=A-C 
Z2=q-n 

IE(A0S(71)  .GT.  0.001  .AND.  AB5(72)  .GT.  0.001)  GO  TO  700 

GO  TO  400  f 

700  CONTINUE  ; 

IE|.AG(N)=0  j 

lREn=lREO*l  ; 

N = M*1 

KMA<  = M _ 1 ! 

lEIIRED  .IT.  LmaX)  go  TO  400  J 

MORE=l  ^ 

DO  500  TRS=1»5  1 

500  RACKSPACE  2 i 

900  backspace  2 'I 

NMAX=1 

K = 1 j 

IE ( IRED  .1  T.  2)  GO  TO  100 
WRITE (6,9000) 

1000  DO  1400  J1=1.Kmax 
DO  1300  11=  1,  2 
ILIKE  = 0 

IE(1ELAG(i1)  .GT.  0)  GO  TO  1400 
A=X(I1,J1) 

R=Y(I1,J1)  ; 

DO  1200  J = 1.  kmax 
DO  1100  1=1,2 

IE(.(  .ED.  (1  .AND.  I .EO.  ID  GO  TO  1200 
1E(IELAG(  I)  .GT.  0)  GO  TO  1200 
AA  = A-X  (1,1) 


KMUiOrfMMaiik 
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PR  = C)-Y(t.  )) 

ir(ABS(Aft)  .GT.  0.001  .ANn.  ARc(RH)  .r,T.  0.001)  GO  TO  1100 
ILIKF  = UTKF  ♦ 1 

GO  TO  noo 

1100  rONTIMUF 
1^00  rONTIMUF 

IFULIKr  .ro.  0)  f,n  To  1700 
noo  CONTlNUr 
1400  co^^Tt^'ur 

no  isoo  ji  = 1.  XMAX 
n = 1 

rr(TrLAG(ll)  .GT.  0)  GO  TO  ISOi 
GO  TO  1600 
1500  r-OMTiKjiir 
1600  A=XfIl,,il) 

R=Y(mji) 

1700  WRITFI7)  FONT. A, R 

WRITE  (6.0)  00)  COX'T.A.R 
K = 1 

IF (T1  .FO.  1 ) K = ? 

WPTTE(7)  rnxJT.X  (X,  Jl ) , Y (K,  ,J1 ) 

WRITE  (6,91  00)  CONT,  X (K  , Jl  ) , Y (K  . ,11  ) 

A = X (K  . Jl  ) 

R = V(K, Jl ) 

TFLAG(J1)=1 

IFMK  - o j 

noo  00  ’000  J=1,kyax 

IFnrLAG(J)  .GT.  0)  GO  TO  ?000 

no  1900  1=1,2  I 

71  =A-X  ( I , ,)) 

7^  = a-Y  (1,1) 

IF(A8S(71)  ,GT.  C.OOl  .OP.  AP';(7?)  ,GT.  0.001  ) GO  TO  1900 
IFLAG(J)=1 
IfHX=0 
N = ? 

IF(T  .'■O.  ’)  >')=1 
A=X J) 

■R  = Y (N,  J) 

WPITF(7)  rOMT.A.H 
WPITF (6.91 00)  rOMT.A.M 
NMA  X = NMA  X ♦ 1 
GO  TO  2000 

noo  coMTiNijF 

2000  CONTINUF 

IF(K|MAX  .rr).  kmaX)  GO  TO  2100 
(►(TCHK  .r,T,  0)  NMAX  = MMAX»1 
nnCHK  .GT.  0)  GO  TO  2100 

1GMK=1 

GO  TO  noo 

2100  COflTlMUr 

r,OMT  = -l . *r0NT 

irdSTOO  .GT.  01  r,o  TO  2^00 

TFnjMAT  K’-tAX  .ANII.  MOPF  .CO.  0)  GO  TO  100 

IF(N)'A'  .to.  Kv-AX  .ANn.  MOPC  .CO.  1)  GO  TO  200 
2100  GO  TO  lOOo 
’400  C0>‘T1N(IC 
PC  Ti  iPN 
’GOO  TGT'^P=1 

GO  TO  1000 

9000  F "ih'^A  T ( T-i  n CorJTOi  |P«..  X » « « 1 i,xo  Y»// ' 

0100  FOP‘-'AT(1h  c 1 2 .4 , 2 ( ?X  .c  1 4. 7)  ) 

F',' 

5S 

J 
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0V‘^OLAV(CnMT0tlP.'?»0» 

PROr-pRAM  moplot 
COmmON/CORF/KOPP 

CORMON/COMTRL/NrOMTRfVAL  (?f»)  . I FLOW , I NUM  t JM  A X »N'ELM,  ! PTS 
COPMON/MXMN/XMAX.XMIN. YMAX, YMiM.xr.YC 
COMMON  A ( n 

CALI  KOPEFL (A.KORF.KFL » 

K0RF  = 6A*  (K0PF/f>4) 
no  100  I=l.KORF 
100  A(T)=0.0 

CAl  L CONNFC  <61 
REWTNn  7 

RFAO(7)  JMAXfWn  M 

LK=  )MAX*1 

LK1=LK*JMAX 

LniF=KORE-LKl-JMAX 

IFd.OIF  .r,T.  0)  00  TO  ISO 

CALL  TFK(l) 

WRITF  (6,<»n00)  LniF,KORF 
CALL  EXIT 
1^0  CONTINUE 

CALL  MXMNPL  ( A (1)  , A (LK)  , A (L*<1)  » IM  A X , SI  7 , IFLOW ) 

DO  ?00  N=1 ,NELM 
,•’0  0 READ  (7)  n ,1^,13 

Mxniim  = 100  This  controls  the  maximum  number  OF  CONTOUR  POINTS 
THAT  will  be  SPLTNED  FITTED  AT  ANY  ONE  TIME 
THIS  MAY  BE  INCREASED  DEPENDING  ON  COPE  AVATl  ABLE 

MXNIlMrlOO 

IK=MXNUM*IPTS 

L<=MXNUM*1 

LKl  =?*MXNllM»l 

LK?=IK*LK1 

LK3=IK+LKP 

LK4=IK*LKT 

LKS=IK*LK4 

LK6-IK*LKS 

LK7=IK*LKS 

LKB=IK*LK7 

LKRsIK+LNO 

LDIF  = K0PE-I.KR 

IFU.DIF  ,OT.  0)  OO  TO  300 

CALL  TEK(l) 

WRITE (FmROOO)  LDIF,K0RE 
CAl  L EXIT 

30  0 CALL  SPLINE  (AM  ) . A (LK ) , A ( LK 1 ) , A ( LK? ) , A (L<3 ) . A (LK4)  ,A  (LKS)  . 

1 A(LKF>).A(IK7).A(LK0)  .MXNUM.SI7) 

9000  FOPMATMH  * CORE  INSUFFICIENT  TN  OVFRl.AY  3«/*  NEEDED  =*I7 
1 * AVAILAPLE  =*I7) 

END 
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SUPOOUTfNr  MXMKJPL  (X.Y,7«JI-'flXtST7.in  OW) 
COMMON/MXMN/XMftX.XMIN.YMAXfYMlKi.xt.YC 
COMMON/HIPLOT/IHI  , ITEt"' 

DIMENSION  X(n,Y(l).7(l>tAl(?).A?(?),lAP(ft) 
DATA  A1  , A?/<.*0.0/ 

X«AX=-1 .E?0 

YMAX=-1 ,E?n 

XMIN=  l.E?0 

Ymin=  1.E?0 

DO  100  J=1,JMAX 

PEAD(7)  X ( J) ,Y ( J) .7( J) 

XMAX=AMAX1 (XMAX.X (J) ) 

XMIig=AMIN)  (XMIN,X  ( J)  ) 

YMAY=AMAX1 (YMAX,Y(J) ) 

YMIN=AMIN1 ( YMIN, Y ( J) ) 

ion  roMTiNijp 

XMAX=AINT(XMAX*0.S) 

XM1N=A1NT (XMlN-O.S) 

YMAX=AIMT ( YMAX+O.S) 

YMIN=AINT<YMIM-0.S) 
lEdTPK  .ME.  1)  GO  TO  TOO 
CALL  TEKd) 

CALL  BINITT 
» CALL  MPTS(?) 

CALL  SLIMX  (POO.BSO) 

CALL  SLTMY(S0,700) 

CALL  OLIMX (XMIN.XmaX) 

CALL  DLIMY(YMIM.YMAX) 

CAIL  XERM(4) 

CALI  YFRM(A) 

CALL  CHECK(Al,A?) 

CALL  0SPLAY(A1 ,A?) 

CALL  CHPSI7(A) 

lEdELOW  .LT.  O)  go  to  ?B0 

00  ?00  J=1 , JMAX 

CALL  FFORm(?(J)  ,(S,1,IAD»32) 

CALL  JUSTER(6.I AR.l . 3? »LEN . 1 OFSFT ) 

CALL  MOVEA(X(J) .Y(J) ) 

CALL  HLABFL(1«A7) 

CALL  M0VRFL(I0FSFT.-5) 

CALL  HLA8FL(LFN.IaR(6-LFN*1) ) 

^00  CONTINUE 
?sn  CALL  BELL 
CALL  TSEND 
CALL  TIMPUT(ICR) 

TOO  IFdHI  .NF.  1)  RETURN 
DX=(XMAX-XMIN)*0.1 
0Y= (YMAX-YMIN»*0.1 

XX=0.0 

YY=n.0 

ANIIM=XMIN 

CALL  NLImBFR(XX.-0.PB«0.  ■ .ANUM. 0.0.1) 

CALL  PLOT (XX. YY. 3) 
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DO  110  1=1.10 
XX=?.*rL04T(I) 

ANUM=ANllM*DX 
CALL  PL0T(XX.YY.2> 

CALL  PLOT (XX. “0.1 ,2) 

CALL  NUMBPR (XX.-0,?5*0.14. ANUM.0.0. 1 ) 
CALL  PLOT (XX.YY.1) 

310  COMTINUF 
XX=0.0 
ANUM=YMIN 

CALL  NUMBPR(“n.?5,YY»0.14.ANljM,q0.«l ) 
CAI  L PLOT(XX.YY,3) 

DO  320  T=1.10 
YY=2.*FL0AT(I> 

ANUM=ANUM+DY 
CALL  PLOT (XX. YY, 2) 

CALL  PLOT(-0.1.YY.?) 

CALL  NIJMBFP  (-0.2S.YY.0.1A.  ANIIM.qO.  .1 ) 
CALL  PL0T(XX.YY,3) 

3?0  CONTIMUF 

517=20. /(ymaX-YmIm) 

IFdFLOW  .LT.  0)  RFTURN 

DO  AOO  J=1 , JMAX 

XX=X(J)*STZ 

YY=Y(J)*ST7 

CALL  PLOT (XX. YV, 3) 

CALL  MAPKFR(l) 

CALL  NIJMBFR(XX*0.1.YY,0.1A,7(J)  .0.0.1) 


400  COMTINUF 


SUBROUTTNF  SPL I NF  ( X * Y , XP , YP . A . p . C . 0 . T . S .MX  , <? I 7 ) 

COMMON/CONTRL/NCOMTP.VAL (?0) , IFLOW.INUM, JMAX.NEI.M, IPTS 
COMMON/HlPLOT/IHl ,ITFK 
COMMON/MXMN/XMAX.XMIN,vmAX. YMIn.XC. YC 

PIMFNSION  X(l).Yn).YP(l),XP(l),A(l),R(l),C(l).0(l),T(l),<=(l) 
dimension  TAR(R) 
so  MORF=0 
IPRMT=0 
IPRNTD=0 
VAK=0.0 
DO  150  1=1, MX 
RFAn(7)  V.Xl.Vl 
TF(F0F(7))  200.100 
100  1F(1  ,FO.  1)  VAK=V 
iio  1F(VAK  ,NF.  V)  SO  TO  190 
X ( I)  =X1 
Y(I)=Y1 
ISO  CONTINUE 
MORF=l 

lon  BACKSPACE  7 
200  1=1-1 

IF<T  .LF,  0)  00  TO  400 
IF (ITEK  ,MF.  I > 00  TO  210 
CALL  FFORm(ABS(VAK) ,6,1 .IAR.32) 

CALL  CHPS17(4) 

CALL  JUSTsR(6,IAR,0,3?,LEN,10FsfT) 

210  CONTINUE 

N0UT=IPTS*I 

J=l 

CALL  CURVE (I ,X,Y,NOUT,XP,YP,A,P,C,n,T,S) 

.Fa,  i>  call  moveaixph > .yph > > 

IFdHI  ,GT.  0)  call  PL0T(XP(l)«SIZ,YP(i)*SI7,3)  j 

21S  J=J»1 

xx=xP( j)*siz 

YY=YP( J) *SI7 

IPRNT=IPRMT*1  J 

IF (MOOllPRNT.INUM) ) 290,220,290  | 

220  IF (ITEK  .FO,  1)  CALL  MOVRFL ( lOFSFT , -S ) I 

IPRNT0=1  I 

IFdHI  ,GT.  0) 

1 call  NUMrfR(XX,YY,0. 21, ABS(VAk), 0.0,1) 

XX1=XP( J-1 )*SI7 
YY1=YP(J-J )*STZ 

IFdHI  .GT.  0)  CALL  PLOT  ( XX 1 , YY 1 , 3 ) 

IFdTEK  .FO.  1)  CALL  HLABFL  (LEm,  I AR  (6-LEN,  1 ) ) < 

IFdTEK  .FO.  1)  CALL  MOVE  A ( XP  ( I- 1 ) , YP  ( J-1  ) ) ] 

IFdPRNTD  .F(3,  2)  GO  TO  50 

POO  IFdTEK  .FO.  1)  CALL  DR  AWA  < XP  ( I)  , YP  ( J)  ) ^ 

IFdHI  .GT,  0)  CALL  PLOT  ( XX  , YY  . 2 ) 

IF(J  .LT.  NOUT)  GO  TO  215 
IF (MORE  .OT.  0)  GO  TO  50 
IFdPRNTO  .GT,  0)  GO  TO  50 

IPPNTD=2  ' 

K=I 

lEdTEK  .FO.  ))  CALL  M0VFA(XP(K)  ,YP(K)  ) 

IFdHI  .GT.  0)  CALL  PLOT(XP(K)*STZ.YP(K)*SI7,3) 

GO  TO  220 
400  CONTINUE 

IFdTEK  ,NF.  1)  return 
CALL  BELL 
CALL  TSENO 
CAi.L  TlNPllTdCR) 

RETURN 

FNO 
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SUPROUTINF  CURVF (NIJMIN»XIN*YIN,NUMOUT«XOUT,YOUT, A,P,C»n.T,F) 
DIMENSION  XIN(l) ,YIN(1) ,XOI)T(l) ,YOUT(l) 

DIMFNSION  Am,R(l),C(l)tD<l).T(l).S(l) 

S(l)  = o.n 
DO  100  T = ?»  NMMIN 
01  = XIN(I)-XINd-l) 
n?  = YIN(I)-YINd-l) 

100  Sd)  = Sd-l)*S0RT(Dl*ni*D?*02) 

400  CALL  SPLCOR(S*XIN,NUMIN,A.R«C.n,T.?..0,«2.0.) 
ns=S (NUMIN) / (NUMOUT-l ) 

DO  SOO  I = 1 f NUMOUT 
TT  = 0S*d-I) 

SOO  CALL  SPLTPP(S. XIN, NUMIN. T«0. TT, XOIJTd)  ) 

CALL  SPLCOF (S.YItJ,NUMIN,A,R,C.n,T,?,0.  ,?.0.) 

DO  FOO  I = I.NIIMOUT 
TT  = DS*d-l) 

600  CALL  SPLTPP(S. YIN, NUMIN, T, 0. TT, YOUTd)  ) 

RFTURN 


DO  n n rt  n n n 


SUPROUTTNF  SPLCOF(X»Y.N*A«B»C»OtT.Il ,V1 .I2.V2) 

THIS  ROUTINE  COMPUTES  THE  VECTOR  OE  SECOND  DERIVATIVES 
IN  ARRAY  T WHICH  are  NEEDED  TO  SPLINE  INTERPOLATE 
TAHliLAR  DATA  STORED  IN  X AND  Y.  ALLOWED  END  CONDITIONS 
INCLUDE  SoECIEICATION  OF  Y"  OR  Y""  AT  EACH  END. 

WHEN  II  IS  1 OR  2 THEN  VI  IS  CHOSEN  VA|  UE  OF 
Y"  OR  Y"",  respectively.  AT  X=X(1).  12  AND  V2  HAVE  A 
SIMILAR  MEANING  EOR  THE  RIGHT  END  AT  X(N). 

DIMENSION  X(l),Y(l).A(n.R(l).rm,D<l).T(l) 

EORM  THE  TRIOIAGOMAL  SYSTEM  DEFINING  ThE  SPLINE 
COEFFICIENTS. 

N1  = N-1 
DO  10  J = 2.  N1 
HJ  = X ( J) -X ( J-1 ) 

HJI=X(J*1)-X(J) 

HJP  = HJ*HJ1 
A(J)  = HJ/HJP 
R(J)  = 2. 

C(J)  = l.-A(J) 

10  D( J)=E.*( ( Y( J*1 )-Y ( J) ) /HJ1-(Y(  l)-Y(J-l) )/HJ)/HJP 
C FORM  THE  CONATIONS  FOP  THE  END  CONDITIONS. 

IE(I1  .EO.  2)  GO  TO  20 
H2  = X (2)-X(ll 
H(l)  =2. 

C(1  > = 1 . 

0(1)  = 6.*( (Y(2)-Y(l) )/H2-Vl)/H2 
GO  TO  30 
20  R(l)=l. 

C(l)  = A. 

D(1 ) = Vl 

30  IE (12  .EO.  2)  GO  TO  AC 
HM  = X(N)-X(N-1) 

A(N)  = 1. 

R(N)  = 2. 

0(N)  = A.*(V2-(Y(N)-Y(N-1 ) )/HN)/HN 
GO  TO  SO 
AO  A(N)  = 0. 

B(N)  = I. 

D(N)  = V2 

C SOLVE  FOR  T(l).  ...  . T(N)  USIhG  GAUSS  REnuCTION 

CO  DO  60  K = 1.  N1 
D = A(K*1 )/R(K) 

R(K*1)  = R(K*1) -R*C (K) 

CO  D (K*l) =D (K*l ) -P*D  (K) 

T(N)  = D(N)/R(N) 

DO  TO  J = 1.  N1 
K = N-J 

■'0  T(K)  = (D<K)-CCX)*T(K*1  ) l/BIK) 

return 

END 
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o o o o r>  o o 


SUHPOUTINF  SPLTRP(X,Y,N.T»ID,XT,FT) 

THIS  ROUTINE  PRODUCES  4 SPLINE  INTERPOl  4TION  VALUE  ET 
CORPESPONDTNG  to  an  X-POSITION  XT.  WHEN  ID  EQUALS 

o»i»2  OP  3f  Then  ET  Takes  the  value  oe'Y,  y's  y**",  or 
THE  integral  OE  Y*DX  PETWEEN  LIMITS  X(j)  AND  X.  THE 
ORIGINAL  DATA  POINTS  ARE  STORE*^  IN  X(I),  Y(I>.  1 = 1, 

...  ,M.  THE  VECTOR  T CONTAINS  THE  Y""  VALUES  RETURNED 
EROM  SUPROUTINE  SPLCOE. 

DIMENSION  X(l),  Y(l),  T(l) 

ini  = 10*1 

C CHECK  whether  XT  IS  OUTSIDE  THr  ORIGINAL  DATA  RANGE 

IE  (XT  .LT.  X(1 ) ) GO  TO  80 
IE (XT  .GT.  X (N) ) GO  TO  120 

C determine  the  first  data  point  to  the  RIGHT  OE  XT 

DO  PO  J = ?,  N 
IE(X(J)  .GE.  XT)  GO  TO  AO 
30  CONTINUE 
AO  J1  = J-1 

HJ  = X(J)-X(J1) 

HJ6  = 6.*HJ 

GO  TO  (=50,80,70,7?)  , IDI 
GO  ET=(T(J1 )*(X(J)-XT)**3  ♦ T(J)*(XT-X(J1) )**3 

% ♦(6.*(Y(J)-Y(J1))-(T(J)-T(J1) )*HJ**2)*(XT-X  ( J1 ) ) 

« ♦8.*HJ*Y(J1)-T(  )1)*HJ**3)/HJ‘> 

RETURN 

<0  ET=(-.3.*T(J1)*(X(J)  -XT)**?*3.*T(J)*(XT-X(J1))**2 
% ♦8.*(Y(J)  - Y(JD)  - (T(J)  - T(J1)  )*HJ**2)/HJ8 

RETURN 

70  ET=(T(J1)*(X(.J)-XT)  ♦ T(  J)*  (XT-X  ( J1  ) ) ) /HJ 
RETURN 

72  HX  =XT-X(J1) 

ET=.2S*T( J1 )*(HJ*«A-(MJ-HX)**A) *.2S*T ( J)*HX**A 
S ♦ (3.*(Y( J)-Y(Jl) ) ♦,S*(T(J1)-T( J) )*HJ**2)*HX**2 
1 ♦(6.*Y(  ll)*HJ-T(Jl)*HJ»*3)*Hv 

ET=ET/HJ6 

IE(  I ,E().  2)  RETURN 
no  7A  K = 2*  J1 
HK  = X(K)-X(K-1) 

74  ET=ET*0.5* ( Y (K-1 ) *Y (K) ) *HK- (T (K-1 ) *T (K) ) *MK*»3/2A. 
RETURN 

aO  H2  = X(2)-  X(l) 
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o o on 


DFTFRMINE  the  slope  at  the  left  eno  ano  interpolate 

LINEARLY. 

SLODE=(V(?)-Y(n  »/H2  - (2.*T(1)  *T(2)  )*H?/(S. 

GO  TO  (RO. 100. 110. 115) , IDl 
QO  FT=Y(1)  ♦ SLOPE*(XT-X(l ) ) 

RETURN 


100  ET=SLOPE 
RETURN 
110  FT=0. 

RETURN 

ns  HX  = XT-X(1) 

ET=Y ( 1 ) *HX*O.S*SLOPE*hX**2 
RETURN 

determine  the  slope  at  the  PIGmT  end  AND  INTERPOLATE 

linearly. 


120  HN=X(N)  - X(N-l) 

SLOPE  = (Y (N) -Y(M-1 ) ) /HN  ♦ (T(N-1)+2.*T(N) )«HN/S. 
GO  TO  (130.100.110.140).  ini 
ITO  ET=Y(N)  ♦ SLOPE* (XT-X (M) ) 

return 

140  HX=XT-X(N) 

ET=Y(N)*HX  ♦ 0.S*SL OPE*HX**2 
DO  150  J = 2.  N 
HJ=X(J)  - X(J-l) 


ISO  ET  = ET*0.5* (Y( J-1 ) .Y ( J) )*HJ-(T  ( )-l ) * T ( J ) ) *H J*« 3/24. 
WFT))RN 
END 
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